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EIGENVALUES  OF  NOPMAL  WAVES  IN  A  PLANE  WAVEGUIDE  CHANNEL. 


G.  I.  Makarov,  V.  V.  Novikov. 


In  preceding  works  [1-2]  the  qualitative  investigation  of  the 
behavior  of  the  eigenvalues  of  the  transverse  operator  of  the  problem 
of  the  propagation  of  electromagnetic  waves  in  the  flat/plane  and 
spherical  waveguide  channels  was-  carried  out.  In  this  case  it  was 
assumed  that  the  lower  (internal)  wall  of  waveguide  is  ideally 
conducting,  and  on  the  upper  wall  impedance  type  approximate  boundary 
conditions  were  used.  It  is  of  interest  to  explain,  to  what  changes 
in  the  behavior  cf  eigenvalues  the  use  of  strict  boundary  conditions 
on  the  upper  waveguide  will  bring.  In  the  present  work  this  question 
for  the  flat/plane  waveguide  will  be  examined.  Since  while  conducting 


of  qualitative  investigations  it  is  necessary  to  have  sufficiently 


simple  analytical  expression  for  the  surface  impedance 


(or  reflection 


coefficient)  of  upper  wall  water/aqueous-water,  appropriate  by  strip 
wave  with  the  assigned,  generally  speaking,  complex  angle  of 
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incidence,  we  will  be  bounded  to  the  examination  of  the  problem  of 
the  propagation  of  electromagnetic  waves  in  the  three-layered  medium 
and  fox  simplicity  we  will  count  as  in  the  preceding/previous  works, 
that  uhe  lower  semi-infinite  layer  is  ideally  conducting.  Analogous 
problem  was  examined  in  the  series/row  of  works  [3-6],  in  which  was 
conducted  the  investigation  of  the  behavior  of  the  eigenvalues  of 
normal  waves  in  some  particular  situations.  The  purpose  of  this  work 
is  a  more  complete  analysis  of  the  characteristic  equation  for 
eigenvalues  and  the  isolation/liberation  of  the  continuous  part  of 
the  spectr  nn. 

§  1.  Representation  of  the  solution  of  problem  in  the  form  of  the  sum 
of  normal  waves. 

.  In  the  present  paragraph  we  will  briefly  reproduce  the 
construction  of  the  solution  of  the  problem  about  the  field  of 
vertical  electric  dipole  in  the  plane-layered  medium,  which  consists 
of  three  uniform  layers,  two  of  which  are  semi-infinite,  and  let  us 
lead  the  necessary  for  further  investigation  transformations  of 
soiut ion . 

Page  4. 


In 


cylindrical  coordinate  system  ir 


:he 


surfaces  of 
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layers  are  described  by  equations  z=0  and  z=h,  and  dipole  is  located 
at  point  r=0,  z=z0  (Fig.  1).  Lower  medi-m  possesses  infinite 
conductivity  (a„  =  oo),  upper  medium  is  isotropic,  uniform  and  possesses 
relative  dielectric  constant  em  and  conductivity  au,and  layer 
irranged/located  between  them,  is  vacuum.  Us  it  will  interest  the 
solution  of  stationary  problem  with  the  dependence  of  current  in  the 
dipole  on  the  time  of  form-  e-,u<.  Mathematical  problem  is  reduced  to  the 
construction  of  the  solution  of  the  equations  of  Maxwell,  who 
satisfies  the  appropriate  boundary  conditions,  including  at  infinity 
(with  z-»<=) .  For  obtaining  the  solution  of  problem  in  the  form  of  the 
series/row  of  normal  waves  it  would  be  possible  to  immediately  seek 
it  in  the  form  of  eigenfunction  expansion  of  transverse  operator.  But 
since  in  this  case  (in  contrast  to  the  approximate  impedance  setting) 
transverse  operator  together  with  the  discrete/digital  possesses 
continuous  spectrum,  that  complicates  construction  the  solution  in 
this  way,  which  is  expedient  to  construct  solution  standardly, 
utilizing  an  expansion  in  the  spectral  functions  of  radial  operator, 
i.e.,  the  integral  transform  of  Four ier-3essel ,  and  then  to  reduce 
him  to  the  series/row  of  normal  waves  and  the  integral  along  the 
continuous  spectrum  of  transverse  operator. 
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!m 


Re  1] 


Fig.  2. 


Page  5. 


Introducing  vector  potential  A=/4e*'and  utilizing  transformation  of 
Fourier-Bessel  during  the  solution  of  the  equation  of  Helmholtz  for 
the  potential,  it  is  not  difficult  to  obtain  the  following  integral 
representation  for  potential  ((Kz<h); 

Hol)  (Xf>  <?(X) 


A  =  &‘. :f  Wt,, 

4it  ^  y ft* -T\*  A  tw  KUK' 


A(A)  (1) 


where  I 
height , 


-  current  in  the  foundation  of  dipole,  hx~  its  effective 


_ QW= _ 

\i  V k%„  ~  Aa  sin  Vk-  -\2(h  —  zu)  — 

-  £m  ]/F^-T2  cos  1/  W  (A  —  z0)l  cos  V" A-  —  Klz  npw  0<z<zo, 
{<  V kh’m  —  Xs  sin  /A8  —  x2  (A  —  2)  — 

[  -  si,,  V A-  —  Aucos  yV  —  X*(A  —  z)]  cos  V  Aa  —  a-z0  rrfw  z0<z<A, 
A  (X)  =  /«*  sin  1/A--XV/  -  j/  A2£m  —  cos  /A--Fa, 


Key:  (1) 
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k  -  wave  number  in  the  vacuum,  =  £,n  +  ^  ~  —  the  complex  relative 
dielectric  constant  of  upper  medium.  In  th<.  expressions  for  Q(X)  and 
A(X)  the  branch  of  funct'-n  Vkh''~^T?  is  f  ixsd/recorded  by  che 

'  rrx 

condition 

1m  ^  —  t?  >  0.  (2) 

As  far  as  function  /k  '-X 2  is  concerned,  fixation  of  its  branch  is 
arbifary,  since  the  integrand  of  integral  (1)  does  not  have  a 
branc'ng  at  prints  X=  k.  The  components  of  electromagnetic  field  are 
expri  .--d  as  po-1  tnt ; al  A  with  the  help  of  the  differential 
relationships/. Ta  ;•  ,  (with  O^r^h) 

dA  F _ E  'i. 

—  dr  ’  '  <u,«o  ozdr  ’  *  <<u*u  r  or\  or  I 

During  the  representation  of  solution  (1)  in  the  form  of  the 
series/r^w  of  .normal  waves  expediently  for  convenience  in  the 
comparison  with  the  results  of  the  problem  in  the  approximate 
impedance  setting  of  passing  in  integral  (1)  to  the  new  variable  of 
the  integration 

7]  = /&-  —  *-A,  Im  V  —  1-  A  =  Irn  0, 

'i 

as  a  result  of  what  (1)  accepts  the  form 


where 
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Q(>  <)  = 

*V sin  7'  ( 1  ~  >f)  ~  V'i ^  Yi ( 1  -  t) 

*V  sim^l  —  jj  —  s^cos^l  —  jj-J  co $*]j, 


COST,  — ,  0<r<z0. 


A  (v|)  ■=  /e,„  •/(  sin  t)  —  kj*  cos  -»),  7j»  =  ]/  f,'-  (sm  —  1 ) ,  $  =  /th, 

in  this  case  Iintj’^o  in  view  of  condition  (2),  and  the  branch  of 

function  j/"  —  pis  fixed/recorded  by  condition  Im  ]/a2  —  di  >  o. 


Page  6. 


The  contour  of  integration  L  goes  along  the  coasts  of  section/cut 
_ _  * 

Im  j/~/r  —  =  o  an<i  AN-  Fig.  2,  where  it  is  given  for  the  case,  is 

shown,  when  wave  number  k  possesses  a  small  imaginary  part  (lm  k>0). 
The  integrand  of  integral  (3)  in  the  lower  half-plane  of  complex 
plane  in),  besides  branch  point  rj=-kh,  has  pole,  caused  by  bullets 
& (i:\. ,  and  it  possesses  one  additional  branch  point  1>  =  V  Kn  ~  1  of 

function  V,  for  unif ormizat ion  of  which  the  sect  ion/ cut  along  line  Im 
7]*=0  is  carried  out  and  branch  Im  77*>0  is  selected. 


The  duct/contour  of 
axis  is  consistent.  In  t 
poles,  which  lie  at  the 
Riemann  surface  for  the 
condition  Im  tj*>0,  and  w 


integration  L  ir  integral  (3)  with  the  real 
his  case  they  will  be  isolated  deductions  in 
lower  half-plane  on  the  upper  sheet  of 
function  rj’f  which  is  fixed/recorded  by 
ill  arise  integral  on  the  duct/contour  G 
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(Fig.  2),  which  encompasses  section/cut  Im  r}*= 0.  Integral  along  the 
real  axis  is  equal  to  zero  in  view  of  the  odd  parity  of  integrand  in 
integral  (3),  so  that  as  a  result  we  come  to  the  following 
representation  of  solution: 


a_13l 


in  which  /„  (z)  =  cos  .  /„  (z„)  =-- cos  ^  —  the  so-called  transverse  functions 

or  high-altitude  factors,  A',-“  excitation  coefficients 

\  — _ ! _  — _ 1 _  /i]\ 

*n  ' 

^  we  the  roots  of  equation  A(TD=U 

"  "  ~  U'nO\  tg  1  =  7,*.  7,*  =  V 7,*  +  ?*  (•;  —  1 )  .  (6)  • 

located  on  the  upper  sheet  of  riemann  surface  for  t»  * 

Expression  (4)  is  the  decomposition  of  solution  in  the 
series/row  by  the  eigenfunctions  of  the  transverse  operator  of  the 
problems,  which  on  the  interval  0Sz<h  coincide  with  transverse 
functions  M2).  and  integral  in  the  functions  of  the  continuous 
spectrum  of  this  operator.  Follows  the  terminology,  introduced  to  ?. 
Ye.  Krasnushkin  [7]  and  L.  M.  Brekhovskikh  [4],  we  will  call 
series/row  according  to  the  eigenfunctions  the  series  of  normal 
waves,  and  integral  term  in  (4)  -  bv  side  wave. 
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Equation  (6)  is  characteristic  equation  for  the  eigenvalues  of  the 
transverse  operator  of  problem,  u,  connected  with  t\  relation  [1] 

t)  =  ,3 1/7- 


In  this  work  we  will  conduct  qualitative  research  of  the 
behavior  of  eigenvalues  in  the  dependence  on  the  parameters  of 
waveguide  and  frequency  of  electromagnetic  vibrations  and  let  us 
compare  results  with  the  results  of  problem  in  the  approximate 
impedance  setting.  In  the  latter  case  transverse  operator  possesses 
purely  discrete  spectrum  [1]  and  expression  for  the  potential  takes 
the  form 

■W=-sr  V]  ]/*  1  -i) .  (7) 

n  o 


where  /-, (z)  and  /.,(-»)  has  the  same  expression,  that  also  in  (4), 


» "‘i’) _ • 

v  n  - 


,  'll'  -‘In 
1 


(8) 


'•  +  < 


'%  are  the  roots  of  the  eauation 


t 


■itti  o  = 


V  *7+ 1 


(9) 

(10) 


Characteristic  equation  upon  a  strict  setting  of  problem  (6)  .for 
the  form  is  reduced  to  the  characteristic  equation  upon  impedance 
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setting  (9),  but  in  this  case  it  is  necessary  to  consider  the 
dependence  of  the  given  surface  impedance  on  the  eigenvalue,  i.e., 
from  the  angle  of  incidence  in  the  normal  wave  on  the  wall  of 
waveguide,  and  to  assume/set 


+ 


(11) 


Approximation  for  surface  impedance  (10)  is  close  to  exact 
expression  (11)  with  satisfaction  of  the  conditions 

«•  ki  »i.  (i2) 

The  investigation  of  the  first  condition  requires  the  knowledge  of 
eigenvalues  u;  we  will  lead  it  after  approximations  for  them  will  be 
obtained. 


Page  8 


In  the  approximate  impedance  setting  the  side  wave,  determined 
by  the  continuous  spectrum  of  transverse  operator  under  strict 
boundary  conditions,  is  absent.  There  should  be  noted  two  limiting 


cases :  s"' —  1  and  h-*0  (with  z=zo=0),  when  the  solution  of  Drobiem  in  a 
strict  setting  is  determined  by  the  exclusively  continuous  part  of 
the  spectrum,  i.e.,  by  integral  term  in  expression  (4).  It  is 
possible  to  show  that  when  \  we  from  (4)  obtain  the  obvious  result 
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R  =  Vr>  +  (z-  zoy- ,  R,  =  v?-  +  (z  +  zuf , 


and  with  h-*0  and  2=2 o=0 


A 1  —  1,1  *  ■  «-  *  ^  ' 

H  1/1-0  —  ~j-~  £,„  - - - 


This  speaks  that  the  impedance  approximate  setting  is  not  valid  at 
the  low  altitudes  of  waveguide  channel,  i.e.,  first  condition  (12) 
must  give  the  specific  limitations  to  the  height/altitude  it  is  deep. 
The  case  of  thin  waveguides  in  detail  is  examined  in  work  [8]. 


S  2.  Investigation  of  characteristic  equation  in  a  strict  setting 


Let  us  consider  characteristic  equation  (6)  on  bipinnate  riemann 
surface  for  function  '(*  =  Kr '-+■  l)  (on  the  upper  sheet  Im  r?*> 0,  on 

the  lower  -•  Im  i?*<0).  The  qualitative  investigation  of  the  behavior 
of  the  roots  of  characteristic  equation  in  the  c*  ipendence  on  and  j3 
in  the  general  case  is  sufficiently  complicated  aJ.gebraic  problem. 

For  simplification  in  the  investigation  we  will  introduce  the 
specific  limitations  to. the  value  of  the  complex  relative  dielectric 


constant 


Sm=£m-Ma,  2==i>0. 

'“tu 


First,  we  will  consider  it  in  the  majority  of  the  cases  carried  out 
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the  inequality 

Ki  ✓  '• 

in  the  second  place,  let  us  asstime 

e*  >  0. 


The  roots  of  characteristic  equation  (6)  on  this  sheet  of 
riemann  surface  are  arranged/located  symmetrically  relative  to 
beginning  of  coordinates;  therefore  it  suffices  to  consider  equation 
in  the  lower  half-plane,  where  Im  rj<0,  especially  because  the 
transformation  of  the  duct/contour  of  integration  in  the  solution  of 
problem  (3)  precisely  only  in  this  half-plane  occurs. 

Page  9. 

Initial  equation  (6)  is  reduced  to  two  real  equations 

a  sill  2a  —  t  ah  -t  _  im  1  in  i^-aRei]* 

2  (ah2 1  cos2  a)  £i  .  ,2  *  U'->) 

m 

sin  2a  1 1 

sx  2  a  -|~  2t  _  _  Kl-  3  I  in  t* 

sll-  t  -j-  CUa-  a  i-n  -j-  a i 

where  a=  Re  tj  ,  r=  Im  rj .  Fig.  3  depicts  the  lower  half-plane  of 

'j 

complex  variable  tj  with  section/cut  Im  rj*=Q  and  signs  Re  r;  and  Im  77, 
which  correspond  to  the  upper  sheet  of  riemann  surface,  are 
indicated.  On  the  lower  sheet  real  and  imaginary  parts  tj  have  signs, 
reverse/ inverse  indicated  in  the  figure. 
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From  equations  (13)  and  (14)  and  that  fact  that  on  the  real  and 
imaginary  axes  signs  Roit*  and  Im  r?  are  identical,  it  follows  that 
with  the  finite  values  /5  and  e,„  the  roots  of  equation  (6)  do  not 
intersect  real  and  imaginary  axes.  Further,  from  equation  (14)  taking 
into  account  the  inequality 


follows  that  the  roots  of  equation  (6)  can  be  arranged/located  only 
the  following  regions  of  lower  half-plane  ( 77 ) : 

1)  in  the  fourth  quadrant  on  the  upper  sheet  of  riemann  surface; 

2)  in  the  fourth  quadrant  on  the  lower  sheet  of  lower  than  the 
section/cut  and  its  continuation; 

3)  in  the  third  quadrant  on  the  lower  sheet  of  riemann  surface. 


Since  in  the  t 
sheet  and  with  the 
not  of  interest  for 


hird  quadrant  the  roots  are  only  on  the 

change  ;3  and  remain  in  this  quadrant, 
'2 

us  and  further  we  can  examine  equation 


lower 
they  are 
(6)  only 


in  the  fourth  quadrant  of  plane  (77).  Let  us  note  the  two  additional 
general/common/total  properties  of  the  roots  of  equations  (6),  which 
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will  be  useful  in  the  examination  of  their  behavior,  roots  can 
intersect  the  upper  shore  of  section/cut  and  the  continuation  of 
saction/cut  on  the  upper  sheet  of  riemann  surface  and  do  not 
intersect  '•lie  lower  shore  of  section/cut  on  the  upper  sheet  and  the 
continuation  of  section/cut  on  the  lower. 
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-fig  3.  Fig.  4. 


Page  10, 

Further  we  will  examine  the  motion  of  the  roots  of 
characteristic  equation  with  that  fixed/recorded  and  change  (3  from 
0  to  «  (i.e.  actually  so  on  change  in  altitude  of  waveguide  h).  With 
[3=0  'sect ion/cut  Im  tj*=0  occurs  along  the  real  axis,  and 
characteristic  equation  has  the  following  roots:  dual  root  in  the 
beginning  of  coordinates  rjo^0  and  simple  roots  on  the  lower  sheet  of 
riemann  surface 

7,^1  =  rnr.  —  arctg  — r—  —  mit — V,  rn—  1;  2,  3,  ...  (13) 

ert»  lm 

'X 

on  che  upper  sheet  of  riemann  surface  roots  are  absent.  The 
qualitative  picture  of  the  location  of  roots  with  (3=0  (in  our  old 
terminology  [1]  -  points,  "output”)  is  represented  in  Fig.  4. 
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With  the  increase  in  /3  with  fixed  value  the 
displacement/movement  of  section/cut  and  roots  of  characteristic 
equation  occurs.  The  equation  of  section/cut  Im  tj*=0  in  the  explicit 
form  takes  the  form 


(16) 


> 


n*  +  «*-(««- 1) 
*) 


whence  it  is  directly  evident  that  for  a=const  to  large  /32  correspond 
the  high  in  the  modulus /module  values  r.  In  this  case  branch  point 
y|b  =  —  if)  yV(  —  |  (Iin  V —  l  0)with  an  increase  £  is  driven  out  from  the 
origin  of  coordinates  along  straight  linearg*i  =  _  +  i.arg(s;(i  _  with 

/3-»®  the  section  goes  away  to  infinity  in  the  lower  half-plane  (we 
count  a>0).  In  this  case  the  roots  of  characteristic  equation,  except 
one,  approach  real  axis  and  have  the  limiting  values 


According  to  the  general/common/total  properties  of  roots  noted  in 

the  beginning  of  paragraph,  these  limiting  values  (point  of  "input") 

lie  on  the  upper  sheet  of  riemann  surface.  3esides  the  entrance 

'2 

points  located  on  the  real  axis,  equation  (6)  has  one  additional 
root,  which  exits  to  infinity  in  the  lower  half-plane  with  /3-*00 


(17) 
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Moreover  it  is  not  difficult  to  show,  utilizing  (16)  and  (17),  that 
it  lies/rests  above  section/cut  and  its  continuation,  i.e.,  it  is 
arranged/located  also  on  the  upper  sheet  of  riemann  surface. 


Page  11. 


Thus,  if  with  ,6=0  all  roots  are  located  on  the  lower  sheet,  besides 
one,  which  lies  in  the  beginning  of  coordinates,  then  with  they 
pass  to  the  upper  sheet,  intersecting  the  upper  shore  of  section/cut 
on  it.  However,  before  analyzing  a  question  about  the  location  of  the 
roots  of  characteristic  equation  relative  to  section/cut,  let  us  make 
some  observation  about  the  course  by  a  line  of  zeros,  whose  equation 
takes  the  form 


1  +  llj*  ® 

1m  T~  =  0. 

1  ”  «m 
1  I L'1 1) 

Re  v  -  >  0. 


since  from  (5)  it  follows 


(18) 


?-=f 


1  TJmJiii-1 

i  — 


The  equation  of  lines  of  zeros  in  case  (18)  in  question  is  much 
more  complicated  chan  upon  the  approximate  impedance  setting  of 
problem  [1];  therefore  here  we  will  be  bounded  only  to  qualitative 
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considerations  and  let  us  give  the  results  of  numerical  calculations 
for  the  specific  cases.  Characteristic  equation  (6)  considerably  is 
simplified  in  two  limiting  cases 


а)  '}  |  *««  “  1  i  <  I  Ti!. 

б)  [i\  -  1  |>hi|. 


In  the  case  of  a)  the  roots  of  equation  are  given  by  expressions 
(15)  and  is  tested  very  small  bias/displacement  with  the 
increase/growth  /3  (see  formula  (32a)).  In  the  case  of  b)  equation  (6) 
for  the  form  is  reduced  to  the  equation  in  the  approximate  impedance 
setting 

=  r  =  -ioy,  (19) 

where  ..  V*'m  —  i 

°*  =  — ~ — •  (2°) 

rrt 

Thus,  in  this  case  (i.e.  the  vicinity  of  entrance  points  lines  of 
zeros  (18)  are  very  close  to  lines  of  zeros  in  the  approximate 
impedance  setting  with  the  surface  impedance  of  the  upper  wall  of 
waveguide,  determined  expression  (20).  If  we  count  then 

v 

expression  (20)  can  be  simplified  and  reduced  to  the  usually  utilized 
expression  for  the  surface  impedance  of  homogeneous  conducting 
semi-infinite  medium  (10) 
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In  this  case  surface  impedance  possesses  the  always  inductive 
character  ( Im  62*<0). 


Page  12. 


However,  relations  (19)  and  (20)  are  valid  in  region  b)  and  when 
inequality  |emj)»ldoes  not  occur.  In  this  case  depending  on 
relationship/ratio  and  a  surface  impedance  6\  can  possess  both 


inductive  and  capacitive  ( Im  §!*>())  character,  namely:  1) 

0  <sm  < 2;  a  <  /•„  (2 - aj  -  Im i]  >  0.  *>  V"6-«(2  —  ««) ~  Ira  o,'<0,  2)  em  >  2  and  a 


arbitrary  (a>0)  -  Im  52*<0. 


Therefore,  generally  speaking,  lines  of  zeros  (18)  in  the  vicinity  of 
the  entrance  points  located  on  the  real  axis,  can  behave  depending  on 
concrete/specific/actual  values  e,„  and  a  as  lines  of  zeros  in  the 
approximate  impedance  setting,  which  correspond  either  to  inductive 
or  capacitive  surface  impedance.  As  far  as  line  of  zeros,  which  exits 
with  the  increase/growth  £  to  infinity  (it  we  we  will  call  infinite 
line  of  zeros),  is  concerned,  it  with«.«>0  and  a>0  always  exisz 
independent  of  the  sign  of  imaginary  part  6*  and  with  the  large  ones 
ji  virtually  coincides  with  lines  of  zeros  in  the  approximate 
impedance  setting  with  the  surface  impedance  6’, determined  by 
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expression  (21).  This  line  of  zeros  is  absent  only  with  a=0,  when 
surface  impedance  6.j  (21)  is  purely  active,  and  surface  impedance  b\ 
proves  to  be  either  also  purely  active  or  purely  capacitive 

(0<t'„,<l).  In  the  case  in  question  is  sufficient  large  ones  j3,  as  we 
see,  surface  impedance  of  the  upper  wall  of  waveguide  proves  to  be 
different  for  the  line  of  zeros,  which  are  terminated  at  the  real 
entrance  points  *0  ,  and  for  line  of  zeros,  which  exits  to  infinity 
(6.*)  .  This  is  caused  by  the  dependence  of  surface  impedance  (11)  on 
the  angle  of  incidence  in  the  waves  on  the  upper  wall  of  waveguide, 
which  must  be  considered,  if  modulus/module  em  is  insufficiently 
great. 


The  study  of  the  behavior  of  lines  of  zeros  in  the  region,  where 
$\V  1 1 ''-hi.  it  is  possible  to  lead  actual  ones  only  by  numerical 

methods;  therefore  we  in  this  case  will  make  only  the  following 
observation.  As  can  be  seen  from  expressions  (15)  and  Fig.  4,  with 
/3=0  the  roots  of  characteristic  equation  are  arranged/located  on  the 
lower  sheet  of  riemann  surface,  except  one,  arranged/located  in  the 
beginning  of  coordinates  and  requiring  separate  examination.  With  the 
increase/growth  (1  the  roots  begin  to  be  moved  along  the  line  of 


zeros,  whose  initial  sections  lie/rest  on  the  lower  sheet  of  lower 


than  the  section/cut  and  its  continuation.  It  is  not  difficult  to  be 


convinced  of  the  fact  that  in  these  sections 


the  surface  impedance 


6(7?)  (11)  has  inductive  character.  In  proportion  to  further  increase 
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/3  the  roots  reach  the  upper  shore  of  section/cut  on  the  upper  sheet 
and  they  emerge  to  the  upper  sheet,  moreover  on  the  upper  shore  of 
section/cut  as  before  Im  5*<0. 

Page  13. 

And  only  when  section/cut  departs  sufficiently  far  from  the  roots, 
which  move  in  turn  to  the  entrance  points  along  lines  of  zeros, 
surface  impedance  5 (77)  can  change  its  character  from  the  inductive  to 
the  capacitive.  Thus,  at  least  in  the  initial  sections  of  the  line  of 
zero  the  imaginary  part  5  ( 77 )  is  always  negative  and  they  are 
arranged/located  in  the  regions,  where  lie/rest  lines  of  zeros  for 
the  inductive  surface  impedance  in  the  case  of  the  approximate 
impedance  formulation  of  the  problem  and  at  which  is  possible  the 
contact  of  lines  of  zeros.  However,  the  contact  of  lines  of  zeros 
leads  .to  the  appearance  of  line  of  zeros,  that  goes  to  infinity, 
which  in  detail  was  considered  in  the  work  of  the  authors  [1]. 

The  phenomenon  of  the  contact  of  lines  of  zeros  is  connected 
with  the  possibility  of  degenerating  the  eigenvalues  of  the 

*1 

transverse  operator  of  problem,  i.e.,  with  the  presence  in  the 
characteristic  equation  of  multiple  root.  In  the  application/appendix 
to  this  work  it  is  shown  that  at  the  real  frequencies, 
examined/considered  by  us,  the  multiplicity  of  the  roots  of 
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characteristic  equation,  i.e.,  the  algebraic  multiplicity  of 
eigenvalues,  cannot  exceed  two,  in  this  case  the  root  of  the  second 
multiplicity  is  determined  by  system  of  equations 

Sill  COS  "»|  Vj  -i.  -  -  ,l"s'  ''  -  (2‘Jti) 

S»*  =  ^  ■  (226) 

For  the  determination  of  the  roots  of  the  second  multiplicity  it  is 
necessary  actual  with  data  s’m  to  solve  equation  (22a)  and  to 
determine  its  roots  (/i  =  0,  1,  2,  ...).  on  obtained  value  7in(£m)  one 
should  then  compute  on  (22b)  the  appropriate  value  of  /32.  The 
degeneration  of  the  eigenvalues  (i.e.  the  root  of  the  second 
multiplicity)  occurs  with  similar  e»'  for  which  the  obtained  value  of 
)32  proves  to  be  positive.  These  values  and  )32  correspond  to  the 
case  of  degeneration. 


Let  us  recall  that  in  the  case  of  the  approximate  impedance 
formulation  of  the  problem  the  multiple  root  satisfies  equation  [1] 

sin  cos  v,-f- *i  =  0,  (23) 

and  value  <  =  —$*»  8  =  -p=====  corresponding  to  degeneration  is  determined 


by  the  expression 


t  =?|tg-»i. 


Pace  14 
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With  fulfilling  of  inequality  1 6„, l  1  the  roots  of  equation  (22a)  are 
close  to  the  roots  of  equation  •»  and  the  approximation  formula 

\  —  ‘C*  ~  ~~~ 5  (uT  ’ 

occurs  moreover  in  this  approximation/approach  relationship/ratio 
(22b)  takes  the  form 


?*=-  +  1 )  Vf  «»  +  — - -<>  tg  (-4)- 


So  that  the  imaginary  part  of  y3 2  (24)  would  become  zero  with  Re  /3 2 >0 , 
the  argument  of  value  +  l  (? ^arg(s'tt -}- 1))  must  satisfy  the 


relationship/ratio 


where 


T  a  —  + 


A"  |  «,«  +  1 1'  ’ 


tg  —  A„e 


In  the  approximate  impedance  setting  in  the  right  side  of  formula 
(25)  .the  second  member  is  absent,  so  that 

C  =  (26> 

and  degeneration  occurs  for  the  specific  arguments  of  value  + 
which  do  not  depend  or  |  e'n  -f-  i  j.  Value  K„  +  1 1  affects  only  the  value  j32, 
which  corresponds  to  degeneration.  As  are  shown  the  results  of  work 

'X 

[1],  value  ’<T"  (26)  for  the  first  numbers  they  prove  to  be  equal  to 

Ko“"  ~  7^"  — 32°,  -t-r" »  24". 

In  strict  statement  when  l!"il  1  the  argument  of  value  +  1  (25), 

which  corresponds  to  degeneration,  depends  onj^-f-li.  In  this  case  the 
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degeneration  at  points  with  n=0  and  n=l  occurs  with  smaller  argument 
fm-f"1  in  comparison  with  the  impedance  approximation/approach,  and  at 
points  n>2  -  with  larger  argument  e'm+  1. 

Let  us  consider  now  the  behavior  of  roots  of  characteristic 
equation  (6) 'and  lines  of  zeros  in  the  vicinity  of  the  entrance 
points  and  output.  In  the  vicinity  of  exit  point  77=0  with 
sufficiently  small  ones  /3,  that  satisfy  the  conditions 


3  |  1/  e‘  —  1  1 

t  j  V  *,«  1  ^  ’ 

(27) 

23  I  ='  V  s'  —  1  v'  1. 

"|  m  '  m  -  ’ 

(2H) 

th*3  root  of  characteristic  equation  is  given  by  the  approximation 

J/V  —  1  (29) 

and  in  this  approximation/approach  it  coincides  with  the  point  of  the 
branching  of  function  y  =  J/V  +  F1*  (em  —  i)  .as  it  follows  from  (29),  in  the 
vicinity  77=0  line  of  zeros,  which  emerges  from  this  point  goes  along 
the  ray/beam,  which  composes  the  angle 

—  ~  T  ~  T  ■  (30) 

with  the  real  axis  a. 

Page  15. 

The  possible  directions  of  line  of  zeros  in  this  case  are  represented 
in  Fig.  5.  When  |£m I  ~  1.  with  not  very  small  ones  /3,  such  that 
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but  simultaneously  and  satisfying  condition  ?  <<  I  it  is  possible  to 
obtain  for  the  root  of  characteristic  equation  in  the  region  1 17 1  <1 
the  following  expression: 


(3t> 


In  this  case  line  of  zeros  is  also  the  section  of  straight  line,  but 


going  with  respect  to  the  real  axis  at  the  angle 


<t>.. 


~  !  'i 

7  -  t  arct«  r-  • 

’  n  fcm 


The  possible  directions  of  line  of  zeros  in  this  case  are  represented 
in  Fig.  6.  Thus,  when  j |  -.>  1  line  of  zeros,  which  emerges  from  the 
point  17=0,  at  first  goes  along  one  of  the  rays/beams  Fig.  5,  and  then 
with  the  increase  in  0  passes  to  the  corresponding  ray/beam  Fig.  6. 
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Fig.  5.  Key:  (1) .  or. 


Fig.  6. 


Page  16. 


In  the  vicinity  of  exit  points  (15)  with  satisfaction  of  the 
conditions 


—  i 


?  C  (wtt) 


for  the  roots  of  characteristic  equation  iz  is  not  difficult  to 
obtain  the  following  approximation: 


1  —  m-  —  arctt'  — 


: '.2 


-(C'+l) 


rnr.  —  arcig  ■ 


(32) 


which  when  K,,!-'5  considerably  is  simplified 
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In  this  case  lines  of  zeros  in  the  vicinity  of  exit  points  (15)  are 
the  line  segments,  which  go  vertically  down. 


With  an  increase  0  the  roots  of  equation  (6)  approach  entrance 
points  rj‘™i  =  |/n  +-rr]- and  with 


for  them  occurs  the  approximation 


(33) 


Lines  of  zeros  in  vicinity  of  exit  points  go  along  the  rays/beams, 
arr anged/located  relative  to  real  axis  a  at  the  angle 


>t> 


+  arctg 


a//  +  s,„«  ’ 


where 


“  +  *'*=  V  e-  -  1  • 

/  L  lf”i  ~  *  jl  -t  +  (fm 


- 1) 


The  course  of  lines  of  zeros  in  the  vicinity  of  entrance  point 

/  .  i  \ 

vji"1'  =  I  m  -f  t)  '  for  different  ones  s,„  and  a  is  shown  m  Fig.  7. 


DOC  =  84100101 


PAGE  27 


The  considerations 

given 

above  and  results  make  it  possible  to  represent  the  dynamics  of  the 
roots  of  characteristic  equation  and  the  behavior  of  lines  of  zeros 
with  different  values  and  a  (£„><).  a>0).  As  a  whole  it  is  possible  to 
say  that  the  general  pattern  of  lines  of  zeros  in  the  case  of  a 
strict  formulation  of  the  problem  in  question  is  the  same  and  in 
impedance  approximation/approach  [1],  with  exception,  separate  parts. 
Corresponding  line  of  zeros  in  the  impedance  approximation  answers 
each  line  of  zeros  in  a  strict  setting. 
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Key  (1) .  or. 

Page  18. 

As  earlier  [1],  it  is  possible  to  number  the  eigenvalues  (roots  of 
characteristic  equation)  on  exit  points  of  lines  of  zeros 
corresponding  to  them.  Lines  of  zeros,  arranged/located  to  the  left 
and  right  of  infinite  line  of  zeros  (17),  go  in  opposite  directions: 
left  lines  of  zeros  go  from  exit  points  into  the  adjacent  entrance 
points  which  lie  to  the  right  of  them,  and  right  lines  of  zeros  - 
into  the  adjacent  entrance  points  which  lie  to  the  left.  Exit  point 
of  infinite  line  of  zeros  depends  on  the  argument  of  value  e,«  and  in 
proportion  to  its  increase  it  is  moved  at  the  beginning  of 
coordinates.  The  displacement  of  exit  point  of  this  line  occurs  at 
the  moment  of  its  contact  with  line  of  zeros,  which  emerges  from 
adjacent  left  exit  point  Numerical  calculations  confirm  the  given 
common  picture  of  the  behavior  of  lines  of  zeros.  As  the  illustration 
Fig.  8  depicts  first  three  lines  of  zeros  for  em=l,5  and  the  different 
values  a.  By  arrovs/pointers  is  indicated  the  direction  of  the  motion 
of  the  eigenvalues  along  lines  of  zeros  with  ar.  increase  /3.  With  cj=0 
(dielectric  waveguide)  infinite  line  of  zeros  is  absent  and  all  lines 
of  zeros  go  in  one  direction  -  from  left  to  right.  With  a=3  infinite 
line  of  zeros  emerges  from  exit  point  of  the  first  mode.  Line  of 
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zeros  for  the  zero  mode  goes  in  this  case  as  before  from  left  to 
right,  and  the  lines  of  all  zero  modes,  beginning  with  the  third,  go 
in  the  opposite  direction. 
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Page  19. 


And  finally  with  a=15  exit  point  of  infinite  line  of  zeros  coincides 
with  exit  point  of  the  line  of  zero  zero  modes,  and  all  lines  of 
zeros,  except  zero,  go  from  right  to  left.  As  was  already  mentioned 
above,  the  displacement  of  exit  point  of  infinite  line  of  zeros  is 
connected  with  the  phenomenon  of  the  contact  of  lines  of  zeros.  This 
phenomenon  for  the  lines  of  zero  zero  and  first  modes  when  «,„  =  1.5  is 
illustrated  by  the  results  of  the  numerical  calculations,  represented 
in  Fig.  S.  Contact  occurs  witha-11.5.  The  argument  of  value  e«  +  1> 
which  corresDonds  to  degeneration  (contact),  comprises 
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approximately/exemplarily  77°45',  which  is  somewhat  less  than  the 
value  of  the  argument  of  this  value  in  impedance 
approximation/approach  (78°). 

Although  the  common  pattern  of  lines  of  zeros  in  a  strict 
setting  is  close  to  rhe  same  in  the  impedance  approximation/approach 
however  on  the  concrete/specific/actual  course  of  lines  of  zeros 
strict  boundary  conditions  can  exert  a  substantial  influence,  which 
most  strongly  is  developed  with  small  ones  |e„J-  Difference  in  the 
concrete/specific/actual  behavior  of  lines  of  zeros  upon  strict  and 
impedance  formulations  of  the  problem  is  connected,  in  particular, 
with  the  following  two  facts,  which  have  already  been  noted  above. 
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First,  exit  points  of  lines  of  zeros,  besides,  zero,  are 
arranged/located  upon  the  strict  setting  in  contrast  to  the  impedance 
not  on  the  real  axis,  their  removal/distance  from  the  real  axis 
depending  on  em  and  a  it  can  be  considerable  (so,  when  and  a-*0 
exit  point  go  away  to  infinity).  In  the  second  place,  in  a  strict 
setting  is  always  an  infinite  line  of  zeros,  which  corresponds  to  the 
inductive  surface  impedance  of  the  upper  wall  of  waveguide  in  the 
impedance  approximation/approach.  However,  in  this  case  in  the  case 
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e,„<-  and  a< V *m(2—  -m)  about  the  entrance  points  of  line  of  zeros  in  a 
strict  setting  they  go  at  sharp  angle  to  the  real  axis  (Fig.  7), 
which  corresponds  in  the  impedance  approximation/approach  to 
capacitive  surface  impedance.  The  latter  fact  leads  to  a  noticeable 
change  in  the  form  of  lines  of  zeros,  especially  lines  of  zeros  with 
the  large  numbers.  Fig.  10  shows  lines  of  zeros  for  e,„=0,J.  and 
different  a,  that  correspond  to  the  ninth  and  tenth  normal  waves. 
Results  given  in  this  figure,  and  also  results  in  Fig.  8  and  9 
sufficiently  clearly  illustrate  the  fundamental  special 
features/peculiarities  of  lines  of  zeros  in  the  case  of  a  strict 
formulation  of  the  problem  in  question. 

Above  we  examined  the  course  of  lines  of  zeros,  without  being 
interested,  on  what  sheet  of  riemann  surface  zero  are 
arranged/located.  Let  us  pause  now  at  this  question,  since  it  has 
vital  importance  in  view  of  the  fact  that  in  the  series/row  of  normal 
waves  for  the  solution  of  problem  in  the  form  (4)  enter  only  the 
terms,  which  correspond  to  the  roots  of  the  characteristic  equation, 
which  are  arranged/located  on  the  upper  sheet  of  riemann  surface  in 
the  lower  half-plane.  Let  us  begin  from  the  examination  of  zero  root, 

'j 

which  corresponds  to  zero  mode  and  which  is  found  with  /3=0  in  the 
beginning  of  coordinates  (r?0^=  0). 
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Fig.  10. 


Page  21. 

With  small  (3,  that  satisfy  inequalities  (27)  and  (28),  this  root 
lies/rests  at  the  vicinity  of  the  point  of  branching  of  function 
( 29) .. Utilizing  equation  (6),  it  is  not  difficult  to  ascertain  that 
it  is  located  on  the  upper  sheet  of  -Siemann  surface,  if  occurs  the 
inequality 

0<*s.  (34) 

also,  on  the  lover  sheet  with  the  reverse/' inverse  inequality 

sm  (e,„  —  1 )  a-’.  (35) 

For  the  zero  root  with  small  ones  £  it  is  possible  to  obtain  the  more 
exact  expression  than  (23),  applying  the  method  of  successive 
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approximations  for  solving  equation  (6),  which  in  the  case  /V<l  in 


question  takes  the  form 


)-<v. 


In  the  first  approximation,  we  obtain 

(*;-!) - 
-  - 1)2- 

From  this  expression  the  following  relationship/ratio  for  the  product 
of  the  real  and  imaginary  parts  ?|(0>  .=  o'°>  + 


23(«)-(«)  =  —  fi-ac  —  4Jj2a  (sm  - 
—  0  — «2|-  (36) 


ensues.  Comparing  relationship/ratio  (36)  with  the  equation  of 
section/cut  and  its  continuation 

.> _ _3*3 

it. is  not  difficult  to  arrive  at  the  following  conclusions  relative 
to  .the  mutual  location  of  root  7j(0'  and  section/cut  or  its  continuation 
(Fig.  11)  in  the  case  of  small  ones  j3: 


when  0  <  e„,  <  —  root  is  arranged/located  on  the  upper  sh.  et.  of 
riemann  surface  in  region  (1c),  limited  by  positive  real  semi-axis, 

’j 

by  the  negative  imaginary  semi-axis,  by  section/cut  and  by  its 
continuation; 


when  — <im<i  root  is  arranged/located  on  the  upper  sheet  in 
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r  gion  (2c),  which  is  arranged/located  in  the  fourth  quadrant  to  the 
right  and  are  lower  than  the  section/cut  its  continuations; 

k,„>1  'and  fulfilling  of  inequality  (34)  root  is  arranged/located 
also  in  region  (2c); 

with  em>l  and  fulfilling  of  inequality  (35)  the  root  is  located 
on  the  lower  sheet  in  region  (2m). 
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In  the  latter  case  with  the  increase/growth  /3  zero  root  passes 
from  the  lower  sheet  of  riemann  surface  to  the  upper,  intersecting 
with  certain  value  3=^  section/cut  Im  7j*=0.  Value  depends  on 
modulus/module  and  argument  of  value Kn =  •+■  a-r ,  <?  =  arctg  ~ J  ■  I n 

order  to  estimate  value  let  us  consider  for  simplicity  the  case, 
when  j  ~'m  |  =p I and  we  will  use  approximation  formula  (31),  valid  with 

Vp  >  ?  »  ■ 

Based  on  this  formula,  it  is  not  difficult  to  show  that  in  case  (35) 
interesting  us  for  the  approximation 


3(0>  = 

•  K|> 


-  *>in  J- 


(37) 


occurs.  The  obtained  exoression  for  is  correct  not  with  very 

•  K|,  1 


ft 


9 
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small1  <p,  such  that 


FOOTNOTE  x.  In  this  case  <p  does  not  exceed  tt/4  (we  count  e"«  ' 1)1  since 
otherwise  zero  root  is  found  on  the  upper  sheet  of  riemann  surface 
already  with  as  small  as  desired  /3.  ENDFOOTNOTE. 


With  the  decrease  of  imaginary  part  by  decrease  ©,  value 
grows/rises  also  with  certain  sufficiently  small  0  expression  (37) 
proves  to  be  wrong.  In  the  case  of  small  ones  t?  the  intersection  of 
the  roots  of  section/cut  can  be  investigated,  relying  on 

1°  I V  l 

approximation  formula  (33),  which  for  zero  root  (m=0)  when  1  1 ' 

takes  the  following  form: 


(38) 


Fulfilling  of  the  inequality 

is  applicability  the  condition  for  this  approximation,  from 
expression  (3S)  it  follows  that  the  root  intersects  section/cut  with 


Formula  (39)  is  valid  for  sufficiently  small  ones  <0 ,  which  satisfy 


the  condition 
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Page  23. 


Thus,  from  expressions  (37)  and  (39)  it  ensues  that  value  l1!';,’*  with 
which  zero  root  transits  from  the  lower  sheet  of  riemann  surface  to 
the  upper,  grows/rises  with  the  decrease  a(-fc— e.  tp)  and  it  approaches 
infinity  with  a-*-0.| 


As  far  as  other  roots  of  the  characteristic  equation  (6), 
arranged/located  in  the  fourth  quadrant,  are  concerned,  they  with  /3= 0 
lie/rest  on  the  lower  sheet  of  Riemann  surface,  and  with  the 
increase/growth  /3  in  the  case  a>0  just  as  zero  root,  they  pass  from 
the  lower  sheet  to  the  upper.  Those  corresponding  to  these  roots  of 
value  with  which  the  roots  intersect  section/cut,  can  be 
determined,  utilizing  approximations  (32a)  (32b)  and  (33).  Thus,  in 
the  case  p>>l  and  sin  <d>>1/t rp  for  %i"’  occurs  the  following 


approximation: 


V  |i“  f  ’ 


from  which  it  foil  that  the  roots  intersect  section/cut  with 
different  ones  whose  value  grows/rises  with  an  increase  in  number 
m.  At  the  same  time  in  the  case  not  of  not  very  small  ones  in 
question  o  zero  root  appears  at  the  upter  sheet  with  the 


.•/ 
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increase/growth  j3  earlier  than  the  first,  which  is  evident  from 
expressions  (37)  and  (40). 

Thus,  with  the  finite  values  of  p  on  the  upper  sheet  of  riemann 
surface  is  located  a  finite  number  of  roots  of  characteristic 
equation  (6),  i.e.,  the  discrete  spectrum  of  the  transverse  operator 
of  problem  consists  of  a  finite  number  of  isolated  points  (it  is 
intended  compare  e,„>0).  With  small  ones  p,  that  satisfy 
relationships/ratios  (27)  and  (28),  the  discrete  spectrum  more  proves 
to  be  empty,  if  the  properties  of  the  upper  medium,  '•  uich  limits 
waveguide,  are  subordinated  to  condition  (35).  Let  us  note  that  in 
the  case  a=0  and  u,>lthe  discrete  spectrum  is  absent  with  any  values 
p.  With  a>0  a  number  of  points  of  discrete  spectrum  depends  on  the 
height/altitude  of  waveguide  (or  value  /3=kh)  and  increases  with  the 
increase/growth  p.  The  fact  of  the  finiteness  of  a  number  of  points 
of  discrete  spectrum  was  not  noticed  in  preceding/previous 
investigations  [3],  [4],  although  for  the  thin  waveguides  it  can  p^ay 
a  significant  role  during  the  representation  of  field  in  the  form  of 
the  series/row  of  normal  waves. 

'2 

§  3.  Some  generalities  about  the  representation  of  the  solution  of 
problem  in  the  form  of  the  series/row  of  normal  waves. 

As  it  was  shown  in  the  preceding/previous  paragraph,  the 
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discrete  spectrum  of  the  transverse  operator  of  problem  when 
consists  of  a  finite  number  of  isolated  points,  i.e.,  in  the 
decomposition  of  solution  in  the  form  (4)  a  finite  number  of  normal 

waves  enters.  For  the  evaluation/estimate  of  number  N  of  normal  waves 

,  ,  .  1 

in  the  case  of  1 j  2>  1  and  iir^  *.n  •'>  _ j 7  j  it  is  possible  to  use 
approximation  formula  (40). 

Page  24. 

When  the  conduction  currents  in  the  upper  medium  considerably  exceed 
bias  currents,  from  formula  (40)  we  find 


where  the  brackets  indicate  whole  part  of  the  number,  included  in 
them1 . 

FOOTNOTE  1 .  Let  us  recall  that  is  the  conductivity  of  the  media, 
which  limits  waveguide  on  top.  ENDFOOTNOTS. 

If  the  height  of  waveguide  is  70  km,  and  the  conductivity  of  upper 
medium  is  equal  to  10' 6  l/ohm*m,  then  an  Earth- ionosphere  in  the 
range  of  the  superlona  waves  corresponds  to  the  effective  properties 
of  waveguide  channel  in  the  daytime  conditions  for  propagation,  then 
a  number  of  normal  waves  at  the  frequency  of  10  kHz  proves  to  be 
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approximately/exemplarily  equal  to  180.  In  this  case  the  fact  that  a 
number  of  normal  waves  is  certain,  does  not  have  vital  importance, 
since  a  number  of  normal  waves,  necessary  for  describing  the  field  at 
the  distances  from  the  source,  large  several  hundred  kilometers,  does 
not  exceed  ten.  So  that  the  number  of  normal  waves  with  the 
properties  of  upper  medium  indicated  would  compose  several  ones,  the 
height/altitude  of  waveguide  channel  must  not  exceed  10  km.  The 
finiteness  of  a  number  of  normal  waves  can  prove  to  be  essential  for 
the  thin  waveguides  or  when  the  medium,  which  limits  waveguide  on 
top,  is  low-loss  dielectric.  In  these  cases  in  the  decomposition  of 
the  solution  of  problem  (4)  the  integral  along  the  continuous 
spectrum  plays  a  noticeable  role. 

The  fact  of  the  finiteness  of  a  number  of  the  eigenvalues-  of  the 
non-self-ad joint  differential  operator  of  the  second  order,  assigned 
on  the  semi-ax.is,  is  sufficiently  familiar  and  investigated  in  the 
series/row  of  mathematical  works  [9]  when  the  coefficients  of 
differential  expression  satisfy  the  specified  conditions.  However, 
obtained  in  these  works  the  results  are  not  directly  applicable  to 
the  examined  by  us  case,  since  there  are  nor  satisfied  the  utilized 

'i 

in  them  conditions.  The  problem  investigated  by  us  shows  that  the 
finiteness  of  a  number  of  eigenvalues  can  occur,  also,  under  the 
wider  assumptions  about  the  behavior  of  coefficients  of  differential 
express  ion . 


alk 
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Expansion  (4)  of  the  solution  of  problem  encompasses  the  series 
of  the  normal  waves  of  those  corresponding  to  the  eigenfunctions  of 
transverse  operator,  and  integral  on  the  duct/contour  G  (Fig.  2), 
which  is  integral  along  the  continuous  spectrum  of  this  operator.  It 
is  here  interesting  to  note  the  following  fact. 

Page  25. 

If  we  examine  the  solution  of  problem  only  in  the  interval  Oizih 
(expansion  (4)  corresponds  to  precisely  this  case),  then  the 
duct/contour  of  the  integration  it  is  possible  to  deform  and  to 
transfer  it  into  the  duct/contour,  whose  part  will  be 
arranged/located  on  the  lower  sheet  of  riemann  surface  for  the 
function  n*(for  example,  duct/contour  G'  in  Fig.  12,  where  by  solid 
line  it  is  shown  the  part  of  the  duct/contour,  arranged,  ocated  on 
the  upper  sheet,  and  by  dotted  line  -  on  the  lower  sheet).  During 
this  deformation  of  the  duct/contour  of  integration  G,  which  can  be 
required,  for  example,  for  computing  the  integral  along  the 
continuous  spectrum  (side  wave)  with  the  help  of  the  method  of 

v 

steepest  descent,  they  will,  generally  speaking,  intersect  the  pole 
of  the  integrands,  which  correspond  to  the  roots  of  characteristic 
equation,  arranged/located  on  the  lower  sheet  of  riemann  surface,  as 
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a  result  in  expansion  (4)  together  with  the  series  of  normal  waves 
will  appear  the  series/row  of  deductions  in  these  poles,  which  we 
will  call  the  series  of  quasi-normal  waves,  since  the  transverse 
functions,  entering  into  them,  are  not  the  eigenfunctions  of 
transverse  operator.  The  transverse  functions  of  quasi-normal  waves 
satisfy  differential  equation  for  the  eigenfunctions  and  boundary 
conditions  with  z=0  and  z=h,  buc  they  do  not  satisfy  boundary 
condition  at  infinity  -  they  exponentially  increase  with  z-*«.  A 
number  of  quasi-normal  waves  in  the  decomposition  of  the  solution  of 
form  (4)  can  be  arbitrary,  including  infinite  (but  calculating), 
depending  on  the  concrete/specific/actual  selection  of  duct/contour 
G'.  During  the  described  deformation  of  duct/contour  G  into  the 
duct/contour  G'  the  field  of  the  side  wave  (integral  along  the 
continuous  spectrum)  is  represented  in  the  form  of  the  series/row  of 
quasi-normal  waves  and  integral  on  the  duct/contour  G',  which 
subsequently  [8]  we  will  call  sky  wave  when  duct/contour  G' 
corresponds  to  the  duct/contour  of  the  fastest  descent,  passing 
through  points  ?,  =  —  i$Y  s'm~~  1  • 

Let  us  pause  now  at  the  properties  of  normal  wave,  which 
corresponds  to  the  eigenvalue,  arranged/located  on  infinite  line  of 
zeros.  The  properties  of  other  normal  waves  are  sufficiently  well 
studied  in  the  series/row  of  works  [3-6].  As  it  was  shown  in  the 
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preceding/previous  paragraph,  among  the  eigenvalues  of  transverse 
operator  with  a>0  is  always  one  value,  which  is  arranged/located  on 
line  zero,  which  exits  to  infinity  with  With  satisfaction  of 

condition  lm this  eigenvalue  is  given  by  approximation 

*  */n 

(17)  . 
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Fig.  12. 


Page  26. 


However,  for  determining  the  excitation  coefficient  of  the 
corresponding  normal  wave  approximation  (17)  proves  to  be 
insufficient  and  it  must  be  made  more  precise.  Utilizing  during  the 
solution  of  equation  (6)  a  method  of  successive  approximations  in  the 
region.  Jim  rjj>l,  it  is  not  difficult  to  obtain  in  zhe  first 
approximation,  the  following  expression  for  the  eigenvalue  in 
question: 


'j/,1 


The  excitation  coefficient  of  this  normal  wave  in  this  case  proves  to 
be  equal  to 
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and  the  part  of  the  vector  potential,  which  corresponds  to  it,  in 
wave  zone  (kr>>l)  whenl^l  '  I  takes  the  form 


V ~sr  ou(r0), 

where  are  introduced  the  following  designations: 


vrm 


+  1 


With  z=z„=h  vector  potential  Alluu  in  this  approximation/approach  is 
given  by  the  expression 

.  .  tlh.,  _ _ 

Aiua  L t,kr~  V *sr  c~'r , 


which  in  the  form  coincides  with  the  ground  wave  of  Tsennek  in  the 
problem  of  Sommerfield  [10]  about  the  propagation  of  the 
electromagnetic  waves  above  the  flat  surface  with  the  source  and  the 
receiver,  arranged/located  on  it.  Value  sr  is  the  numerical  distance, 
introduced  in  this  problem.  High-altitude  factors  of  the 

considered/examined  normal  wave  actually  exponentially  grow/rise  with 
an  increase  in  altitude,  which  indicates  the  fact  that  its  energy  in 
essence  is  concentrated  in  the  vicinity  of  the  upper  wall  of 
waveguide,  in  view  of  these  facts  we  call  this  normal  wave  surface. 


At  conclusion  of  this  paragraph  let  us  consider  applicability 
conditions  for  impedance  of  approximation  during  the  determination  of 
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the  eigenvalues  of  transverse  operator.  As  already  mentioned  in  the 
first  paragraph,  these  conditions  are  determined  by  inequalities 
(12),  the  first  of  which  can  be  rewritten  as  follows,  utilizing  a 
relat ionship/ratic  rt= /fyV: 

M«P2b;-i|-  (42) 

Page  27. 

If  we  utilize  for  t\  approximation  (31),  then  it  is  not  difficult  to 
be  convinced  of  the  fact  that  inequality  (42)  for  the  eigenvalue  of 
zero  mode  in  the  case  |^|>lis  fulfilled  under  the  condition 


(43) 


The  eigenvalues  of  other  normal  waves,  except  one,  corresponding 
to  infinite  line  zero,  possess  the  real  part,  which  exceeds  several 
times  imaginary,  as  is  evident  from  the  results  of  the  numerical 
calculations,  represented  in  Fig.  8-10.  In  this  case  the  real  part  of 
the  eigenvalue  of  m-th  normal  wave  is  of  the  order  m-,  so  that 
inequality  (42)  in  this  case  takes  the  form 


?  »  r^r  ('«=  l.  3>  ■  ■■)• 


(44) 


However,  as  far  as  the  eigenvalue,  which  lies  on  infinite  line  of 
zeros,  is  concerned,  as  it  follows  from  approximation  (17), 
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inequality  (42)  occurs,  if 1 4 1 »  1-  Thus ,  applicability  conditions  for 
the  impedance  approximation/approach  for  eigenvalues  are  determined 
by  inequality  |C|  1  and  inequalities  (43)  and  (44)  respectively  for 

the  zero  and  following  normal  waves.  With  satisfaction  of  these 
conditions  the  excitation  coefficients  of  normal  waves  in  impedance 
approximation/approach  (8)  with  the  high  degree  of  accuracy  coincide 
with  the  excitation  coefficients  in  strict  setting  (5).  Otherwise 
impedance  approximation/approach  leads  to  the  essential  difference  in 
the  properties  of  the  normal  waves  in  comparison  with  the  strict 


setting.  Thus,  for  instance,  for  the  zero  mode  in  the  case  of  the 
thin  waveguide,  when  is  fulfilled  inequality  4(JS  j ««  K,  —  I )  |  l ,  upon  a 


strict  setting  the  eigenvalue  is  determined  by  approximation  formula 
7C  ~  V  —  i ,  whereas  in  the  impedance  approximation/approach  you 


have  yOn 


V  —ii$,  — 


v*«+ 


in  this  case  an  excitation  coefficient  of 


r  ,,  /  ,1  X 

zero  mode  in  a  strict  setting  proves  to  be  equal  to  I). 


and  in  impedance  approximation/approach 


Page  28. 


Inequalities  (43)  and  (44),  are  determined  the  applicability  of 
impedance  approximation/approach,  superimpose  definite  requirements 
for  the  height/altitude  of  waveguide.  For  the  thin  waveguides  the 
impedance  approximate  setting  can  prove  to  be  unsuitable,  in  spite  of 


fulfilling  of  inequality  |=m  »1-  In  the  case  of  waveguides  with  the 
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sufficiently  high  altitude  the  impedance  approximation/approach  is 
inapplicable  for  the  normal  waves  of  high  order  (with  the  high  value 
of  m) ;  however,  these  normal  waves  possess  high  attenuation,  since 
they  prove  to  be  those  not  running,  and  are  unessential  for  the 
description  fields  already  at  comparatively  small  distances  from  the 
source. 

APPENDIX. 


Investigation  of  the  multiplicity  of  the  roots  of  characteristic 
equation. 


Let  us  show  that  characteristic  equation  (6)  does  not  have  roots 
of  the  third  multiplicity  at  the  real  frequencies  (i.e.  with  02>O). 
The.  roots  of  the  second  multiplicity  must  satisfy  the  equations 


eliminating  V  -r  U„  -  M  f  rom  (2n)  with  the  help  of  (Iri)  expressing  /32 
from  (In)  through  and  17,  we  come  to  two  equations  for  the  twofold 
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roots 


CObJ  r. 

bln  <|  COb  /"(  -j-  ij  =  —  — — - 

C  5|»  ’! 

y,  1  ±  1 


t4") 

(on) 


If  we  now  exclude  Vr,{-  - 1)  from  (3n)  with  the  help  of  (In)  and 

with  the  help  of  (4n),  then  equation  (3n)  is  reduced  to  the  following 
form: 

3  bin  >|  cos  >|  +  ( 1  +  -  sin-  <;)  =  o  (on) 

utilizing  this  equation,  it  is  not  difficult  from  equations  (4n)  and 

•  2 

(5n)  to  find  the  expressions  for  em  and)*-..  which  correspond  to  the 
roots  of  the  third  multiplicity 

,2  2  -  cob  2 *  ,  27  <4  +• 1 

-  ~  -"ob  2 1)  -  ?  =  -  T  •  0") 

Thus,  the  roots  of  the  third  multiplicity  are  the  roots  of 
equation  (6n),  which  can  be  rewritten  in  the  following  form: 

/••  (5)  _=  3  bin  i  •+■  ;  (2  ~  cos  ;)  =  0.  -  =■  2<|,  («u) 

but  values  em  and  corresponding  to  them  are  determined  by 
expressions  (7n).  Now  it  is  necessary  to  find  the  roots  of  equation 
( 8 n )  and  corresponding  to  them  values  4  and  /3:.  Equation  (8n)  is 
reduced  to  two  real  equations 

3  bln  x  cli  y  -f  2x  —  x  cob  x  ch  y  —  y  bio  x  bh  y  =  0.  (9n) 

3  cos  -<bh  y  +  ‘2y  —  ycosx  bh  y  +  x  >ln  .c  bh  y  =  0,  (I0u> 


where  ^x-^iy. 
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Page  29. 

Roots  of  equation  (8n)  are  arranged/located  symmetrically  relative  to 
x  and  y  axes;  therefore  it  suffices  to  consider  them  only  in  the 
fourth  quadrant  of  plane  (£). 

Equation  (8n)  has,  in  the  first  place,  a  root  £=0(x=0,  y-0),  for 
which  hml300  and  /3 2  =0 .  This  root  does  not  represent  for  us  interest. 
In  the  second  place,  equation  (8n)  has  the  imaginary  root,  which 
satisfies  equation  (lOn)  with  x=0 

3  ——  =»  ch  y  —  2.  (l  In) 

The  numerical  solution  of  equation  (lln)  leads  to  the  following  value 
of  this  root  and  values  i'„.  and  P2  corresponding  to  it: 

i«  —  — /3,436;  sjo'1  =  ‘>.965;  = -1.530; 

* JiT*  —  —0.965;  ^")2= -0,372. 

(£?n  the  real  axis  (y=0,  x>0)  equation  (8rr),  which  is  led  to  the 

form 

sin  .c 

•J  -j-  3  — - - cos  x  =  0, 

'2 

does  not  have  solutions. 

For  the  localization  of  the  position  of  complex  roots  of 


ft 


equation  (8rt)  and  explanation  of  their  multiplicity  it  is  possible 
to  use  the  Cauchy  theorem  and  a  number  of  roots  of  regular  function 
jjlJ.  Applying  this  theorem  to  function  F(^)  in  the  regions  of 
complex  plane  (§),  limited  by  the  sections  cf  straight  lines  y  =  0, 
y=-C(C»1),  x=C^,  xrC^)  is  not  difficult  to  show  that  the  roots 

of  equation  (8n),  besides  those  noted  above,  are  arranged/located 
in  the  bands 

(2*_T)s<y<5tt*  (£=1,2,3,...), 

moreover  all  roots  are  simple.  During  the  solution  of  system  of 
equations  (9n)— ( 10n)  it  is  expedient  instead  of  x  to  introduce 
a  new  variable 

tk=*2kx—x, 

which  varies  in  an  interval  from  0  to  The  system  of  (9n)-(10rt) 

equation  is  convenient  to  record  in  the  form,  suitable  for  using 
this  method  of  successive  approximations, 


Sln2T 

ch  J'*==2  +  4^ri-_(3chyJ,-y*shy*)-^., 
sin  t*  =*  —  —  —  sins  y  -f-  (3  ch  n  _  ykSh  yk)  sly*  ^ 

jr*  =  2£r.  —  s*.  * 


(!2n) 


Since  **»!  and  2£*;»‘*.  in  the  zero  approximation  we  obtain  the 
the  following  expressions  for  the  roots  of  equations(8rt) : 


4°) « if)  +  /ym  _  ,(0)  f 


ch  y jjj0* 


,(0) - 

*  ~  2 Ax' 

-2,  yW=yP), 


—1,3170. 


(13n) 

(14n) 
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Table  I. 


•  1 

>■* 

‘k 

o 

•* 

(TOMM.) 

U)  R. 

(npHftj.  (15  n',) 

.  ’ <-> 

Jm  iMie  | 

ii>  (TOMH  )  I 

—1,280 

0,5173 

0,4775 

-0,7822 

-0,7990  - 

0,4594 

0,4870 

-1,309 

0,2432 

0,2387 

—0,5066 

-0,5098 

0,3954 

0,3995 

-1.314 

.  0,1605 

0,1592 

-0,4001 

-0,4013 

0,3398 

0,3412 

-1,315 

0,1199 

0,1194 

-0,3403 

-0,3410 

0,3013 

0,3020 

-1,316 

0,09575 

0,09549 

-0,3012 

-0,3015 

0,2733 

0,2736 

-1,316 

0,07972 

0,07958 

-0,2729 

-0,2731 

0,2517 

0,2519 

—1,316 

0,06830 

0.06821 

-0,2512 

-0,2514 

0,2344 

0.2346 

-1,317 

0,05975 

0,05969 

-0,2341 

-0,2342 

0,2203 

0,2205 

-1,317 

0,05310 

0,05305 

-0,2200 

-0,2201 

0,2085 

0,2086 

-1,317 

0,04778 

0,04 « 1 5 

-0,208! 

-0,2082 

0,1983 

0,1984 

—1,317 

0,04343 

0,04341 

-0,1980 

-0,1981 

0,1895 

0,1896 

-1,317 

0,039S1 

0,03979 

-C.1894 

-0,1894 

0.1319 

0,13i2 

-1,317 

0.03674 

0,03673 

-0,1817 

-0,1817 

0,1750 

0,1750 

-1,317 

0,03412 

0,03411 

-0.1748 

-0,1748 

0.1 68? 

0,1689 

-1,317 

0,03184 

0,03183 

-0,1587' 

-0,1687 

0,1633 

0,1633 

-1,317 

0,02985 

0,02954 

-0,163! 

-0,1631 

0,1583 

0,1583 

-1,317 

0,02809 

0.02S09 

-0,1581 

-0,1581 

0,1537 

0,1537 

-1,317 

0,02653 

0,02653 

—0,1535 

—0,1535 

0,1494 

0,1494 

-1,317 

'  0,02513 

0,02513 

-0,1494 

-0,1494 

•  0,1456 

0,1456 

—1,317 

0,02388 

0,02388 

—0,1455 

—0,1455 

0,1420 

0,1420 

(1).  (it  is  precise.).  (2).  (appr.  (15  p). 
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Table  II. 


. ! 

"Sf” 

(l)W*  ■ 

mm 

|7Tgn|imn 

m 

mm 

mtum 

wmm 

m 

Hg 

l 

5,036 

0,5656 

_ 

5,454 

i 

—  17,51 

— 

o 

24,60 

- 

4,604 

— 

45,23 

— 

—  41,90 

-  61,72 

-  89.10 

3 

4 

60,18 

112.8 

1 !  ,58 

21,13 

— 

108.3 

194.3 

121,8 

208,7 

-  68,40 

-  97,31 

5 

6 

7 

182,2 

269,4 

374,8 

— 

32,85 

46,58 

62.27 

~ 

302,0 

431,9 

584,4 

317,7 

448.6 

601,1 

—  128,2 
—161,2 
-196.2 

—118,9 

-150,8 

—184,6 

8 

9 

10 

497,6 

638,2 

798,1 

— 

79,63 

98,59 

119.3 

— 

iOi  ,/ 

958,3 

1170 

775,3 

970,8 

1188 

-232,7 

-270,9 

-310,7 

—220,2 
-257,6 
-  296,5 

11 

975,4 

962,7 

141,2 

126,2 

1407 

1425 

-352,0 

—336,9 

12 

1168 

1 157 

164,3 

148,9 

1662 

1685 

—32  1,0 

— 3 < 3,8 

13 

1382 

1371 

189,2 

172,9 

1943 

1965 

-438,0 

—422,0 

14 

1515 

1602 

215.2 

198,3 

2245 

2267 

-483,4 

— 466.5 

15 

16 

1865 

2138 

1852 

2120 

242,1 

271,4 

225,0 

252,9 

2566 

2913 

2589 

2933 

-530,3 

—577,9 

—512,4 

—559,4 

IT 

2425 

2407 

301,0 

282,0 

3277 

3297 

-625,8 

— 607,6 

18 

2734 

2713 

332,1 

312,3 

3666 

3683 

i  —677,2 

— 65,  ,1 

19 

3037 

363,0 

343,6 

4064 

4088 

-728,1 

— 1 0  <  ,6 

'3398 

3380 

396,0 

376,1 

4492 

i  4515 

1 

t 

-780,3 

t 

-759,2 

1 

Key:  (1).  (it  is  precise.).  (2).  (appr.  p) . 
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In  the  same  approximation/approach  of  value  •«»  and  /3 2 ,  that 
correspond  to  triple  root,  are  given  by  the  formulas 
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./r-jyr.i  /j -»‘°>  1 

y  4/f  it  is  Ij 

fr-)J =» kw  | 

[' +  |/ 

/rj/T  /T- 
y  ikn  is 

.  (I3ii) 
("!.«) 
■  (*"'0 


In  "fable  I  are  given  values^y^  and  obtained  by  the  numerical 
solution  of  system  of  equations  (12n),  and  values  «S,,('3a).  and  also 
precise  values  and  equations  (8n)  calculated  by  approximation 
formula  (15n)  for  the  first  twenty  roots.  As  can  be  seen  from  the 
given  numerical  results,  the  approximation  formulas  for  yk.  £*and  e,’«« 
provide  high  accuracy  already  with  not  very  high  values  of  k  (order 
of  several  units).  The  results  of  the  numerical  calculations  of 
values  p) ►,’<“>*  are  represented  in  Table  II.  The  accuracy  of 
approximation  formulas  (16n)  (17n)  for  K  especially  for  im?l+,*.are 
much  lower  than  formulas  (13n)-(15n)  for  y*.  <,*,  however,  it  proves 
to  be  not  less  than  several  percentages  with  k=20.  from  the  given 
* esults  of  numerical  calculations  for  k=l-20  and  from  approximation 
formulas  (16n)  (17n)  for  k>20  it  follows  that  for  the  roots  of  the 
third  multiplicity  of  initial  characteristic  equation  (6)  value 
has  complex  value.  Consequently,  with  /32>0  (i.e.  real  k=w/c)  the 

multiplicity  of  the  roots  of  characteristic  equation  cannot  exceed 
two. 
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ELECTROMAGNETIC  FIELD  IN  A  PLANE  THTN  WAVEGUIDE. 

G.  I.  Makarov,  V.  V.  Novikov. 

In  the  present  work  electromagnetic  fields,  excited  by  vertical 
electric  dipole  in  the  three-layered  medium  with  plane  interfaces 
(Fig.  1),  will  be  studied.  It  is  assumed  that  lower  half-space  (z<0) 
possesses  infinite  conductivity,  the  properties  of  layer  0Sz£h 
coincide  with  the  properties  of  vacuum,  and  upper  half-space  (z>h) 
represents  homogeneous  medium  with  dielectric  constant  ,  and 
conductivity  a.  Value  can  take  both  positive  and  negative  values 
The  latter/last  case  is  realized,  when  upper  half-space  is  formed  by 
cold  plasma.  In  its  setting  this  problem  coincides  with  the  problem, 
examined  in  the  preceding/previous  work  of  the  authors  [1],  where  th 

general/common/total  properties  of  the  eigenvalues  of  transverse 

“2 

operator  L-  were  investigated. The  problem  about  the  propagation  of 
electromagnetic  waves  in  the  flat/plane  waveguide  was  examined  by 
many  authors,  in  particular  L.  M.  Berkhovskikh  in  monograph  [2]. 
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However,  in  this  work  the  series/row  of  thin  and  interesting,  from 
our  point  of  view,  effects  either  is  not  examined  completely  or  it  is 
examined  inaccurately.  The  fact  is  the  main  result  of  our 
preceding/previous  investigation  [1]  that  the  discrete/digital  part 
of  the  spectrum  of  operator  L.  consists  in  the  case  e,„>0  of  a 
finite  number  1  of  the  eigenvalues,  the  algebraic  multiplicity  of 
each  of  which  can  reach  two. 

FOOTNOTE  x.  Not  denumerable,  as  it  was  assumed  in  works  [2],  [3]. 
ENDFOOTNOTE. 

A  number  of  points  of  discrete  spectrum  (each  of  which  is  compared 
with  the  normal  wave)  depends  on  the  properties  of  upper  half-space 
and  thickness  of  waveguide  h.  It  is  calculating  only,  if  *,,,<0  (this 
case  in  the  work  [2]  it  was  not  examined).  But  if  the  number  of 

normal  waves,  which  participate  in  the  representation  of  fields,  it 
is  always  certainly  and  the  less,  the  thinner  the  waveguide.  With  h 
sufficient  small  ones  the  discrete/digital  part  of  the  spectrum  can 
be  in  all  of  one  eigenvalue  (i.e.  in  the  decomposition  it  is  present 
only  one  normal  wave)  or  even  be  empty. 
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However,  in  the  case  of  the  approximate  impedance  setting  1  the 
spectrum  is  purely  discrete/digital  and  without  depending  on  the 
sizes/dimensions  of  waveguide  and  properties  of  upper  half-space 
consists  of  a  denumberable  number  of  isolated  points. 

FOOTNOTE  1 .  Ia  this  setting  the  dependence  of  the  given  surface 
impedance  5  on  the  eigenvalue  of  transverse  operator  is  ignored. 
ENDFOOTNOTE. 


The  properties  of  normal  waves  for  this  case  are  sufficiently  well 
studied  with  the  help  of  direct  variational  methods  [4].  Thus  with 
small  h  the  impedance  setting  proves  to  be  inapplicable,  and  for  this 
very  reason  the  study  of  thin  waveguides  is  of  special  interest. 
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S  1.  Representation  of  the  solution  in  the  form  of  the  sum  of  normal, 
quasi -normal  and  sky  waves. 


In  order  to  avoid  the  appearance  of  bulky  formulas,  we  will 
subsequently  consider  that  the  source  and  observation  point  are 
arranged/located  on  the  interface  with  lower  half-space  (z=zo=0).  If 
we  use  now  formula  (1)  of  work  [1]  and  to  switch  over  from  the  vector 
potential  A  to  the  components  of  electromagnetic  field,  then  it  will 
seem  that  vertical  component  of  electric  field  takes  the  following 
form:  1 


(l.D 


k  -  wave  number  in  the  vacuum,  and  P„  -  complete  dipole  moment  of 

' 2 

source,  connected  with  current  I,  introduced  in  [1],  by  the 
relationship/ratio 
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FOOTNOTE  2.  We  assume  that  with  the  sufficiently  large  |A|  let  us 
assume  that  the  integration  occurs  in  the  first  quadrant  along  beam 
a=.|X|£/*,  e>o,and  secondly  x  =-|X|e-‘‘.  ENDFOOTNOTE. 

Fixation  of  branch  /k 1 -X 1  is  arbitrary,  on  plane  (X)  from  the  points 
±kj  the  sections/cuts  along  the  lines 

Im]/"  kl-V  =0  (1,3) 

are  carried  out  and  it  is  assumed  that  on  the  upper  sheet 

(1,4) 


Page -35. 

Expansion  (1.1)  contains  both  the  continuous  and  discrete/digital  of 
the  part  of  the  spectrum.  As  in  [1],  we  will  identify  integral  (1.1),  . 
propagated  over  both  shores  of  section/cut,  with  the  expansion  in 

terms  of  the  continuous  part  of  the  spectrum,  and  the  deductions  in 

'2 

zero  denominators  (1.2),  which  are  isolated  with  the  coincidence  of 
duct/contour  (1.1)  with  the  duct/contour,  which  encompasses 
section/cut,  with  the  normal  waves  (or  eigenfunction  expansion  of 
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transverse  operator).  The  case  of  multiple  poles  will  correspond  so 
that  in  expansion  there  participate  also  the  associated  functions  1 . 

FOOTNOTE  1 .  In  work  [1]  it  is  shown  that  with  the  real  u  the 
algebraic  multiplicity  does  not  exceed  two.  ENDFOOTNOTE. 

Together  with  the  section/cut,  let  us  lead  on  both  copies  of  plane 
(A)  from  the  point  A=kj  the  half-lines,  whose  equation  (see  Fig.  2) 
has  the  form: 


Re  X  =  Re  k2, 

(1,5) 

Let  us  assume  that 

IV|>1,  (1.6) 

then  on  the  entire  half-line  (1.5)  it  will  be  possible  to  use  the 
simplest  asymptotic  representation  of  the  Hankel  function  it  will 
seem  that  the  half-lines  (1.5)  present  the  duct /contour  of  the  most 
rapid  descent  for  integral  (1.1).  Therefore  for  the  determination  of 
asymptotic  representations  (1.1)  according  to  parameter  (1.6)  we  had 

v 

to  combine  the  way  of  integration  in  (1.1)  with  the  half-lines  (1.5). 
In  this  case  we  will  cross  certain  number  of  poles  of  expression 
(1.2)  and  will  be  isolated  a  finite  or  denumberable  number  of  the 
terms  outside  the  integral.  One  should  stress  that  obtained  thus 
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representation  of  field  will  no  longer  be  expansion  in  the  form  of 
the  sum  of  normal  waves  and  integral  on  the  continuous  part  of  the 
spectrum,  sin-.e  into  it  the  part  of  the  deductions  in  the  poles, 
arranged/located  on  the  lower  sheet,  will  enter  also,  and  they  are 
not  the  eigenfunctions  of  transverse  operator.  Therefore  we  will  call 
this  sum  the  sum  of  normal  and  quasi-normal  waves,  and  the  component 
which  contains  integral  on  the  duct/contour  of  descent  -  by  sky  wave. 

Page  36. 

$  2.  Investigation  of  the  poles  of  function  F(X,  h) . 

In  this  paragraph  we  will  attempt  to  visualize  the  behavior  of 
zero  equations 

Vki- tg  V  k--\2m  h—-- l—]/  —  .  (2,1) 

m 

depending  on  the  height/altitude  cf  waveguide  h.  The  corresponding 
investigations  on  plane  (rj) 

7)m  =  i/'/c2  -  ,?m  h  ‘  (2,2) 


for  8,n>0  were  made  in  the  preceding/previous  work  of  the  authors 
[1].  In  this  case  us  they  will  interest  mainly  thin  waveguides. 
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Therefore  here  it  suffices  to  study  only  the  rough  picture  of 
behavior  More  detailed  information  will  be  required  only  for  case 
of  h-+0  (i.e.  exit  points  in  our  old  terminology). 

With  h-»0  it  can  be  presented  two  possibilities: 


1 )  |  ]/  k-  —  k 

ih  |  -/-o. 

(2,3) 

A-c  0 

2)  |  Y  k'  -  X 

I  h  |  -  0. 

(2,4) 

h-f  0 

first  case 

i  .  i 

(2,5) 

A-.0 

m  —  0,  ±  1 , 

±  2  ... 

The  location  of  root  xm  on  the  upper  or  lower  sheet  of  riemann 
surface  depends  on  the  sign  of  the  imaginary  part  of  second  term  in 
(2.5).  Let  us  introduce  the  designations 


'  _  /? 

—  p€  i 


(2,6) 


arctg - i=  a  -f-  iv. 


(2,7) 
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VS 


Fig.  2. 


Y=o  (j>*o 

(arr.tg  7^) 

<p-7T /4 

/hK 

0 

7C/2  a 

Fig.  3. 


Page  37. 

If,  furthermore,  to  place  for  the  certainty 


o 


U 


(2,8) 


that  by  comparatively  simple  computations  it  is  possible  to  find  the 
form  of  the  half-lines 

0<p  <  oo,  -f  =  const  (2,9) 

on  complex  plane  (2.7).  The  corresponding  lines  are  represented  in 

'2 

Fig.  3,  the  arrows/pointers  indicating  the  direction  of  the  motion  of 
point  on  plane  (2.7)  with  the  increase/growth  p.  One  should  have  also 
in  mind  that  the  physically  realizable  structures  correspond  to 


DOC  =  84100102 


PAGE 


change  <p  within  the  limits 


0<<p  .'it. 


(2,10) 


It  is  now  not  difficult  to  determine,  in  what  case  Xm  with  h-+0  (exit 
point  is  located  on  the  upper  sheet  of  plane  (X).  For  this  purpose 
let  us  note  that  if  are  satisfied  the  conditions 


rrv c. —  arctg  — r—  >  kh, 

^  m  I 

I  nvx  —  arclg  — I—  >  6  j )/  £,'n  —  1  k , 


(2,11) 


then  \m  it  is  possible  to  find  from  the  expression 


±  -i- |*>  +  i  (/««  —  «)! . 


(2,12) 


Let  us  substitute  now  (2.5)  and  (2.12)  into  the  left  side  of 
transcendental  equation  (2.1),  then  it  seems  that 


Vr-i  * - 


h  •  il  in 


l  me  1  ,  <  1 

-  T  -~arctg—  •  ^ 

in  m 


Consequently  (2.1)  it  will  be  satisfied,  if 


run.  I  ,  ; 

— arctg  ^  . 


(2,14) 
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From  condition  (1.4)  and  formula  (2.7)  it  follows  that  (2.14)  has  the 
place  on  the  upper  sheet  of  plane  (X)  in  the  case  of  v<0,  i.e. 

Reem<0.  (2,10) 


It  is  possible  also  to  be  convinced  of  the  fact  that  with  fl  00 
those  corresponding  to  the  "locked  branches"  [1],  approach  point 
X=±k,  and  corresponding  to  "infinite  branch",  to  the  point 


X  = 


y  4  + 1 


Page  38. 

Since  the  examination  of  thick  waveguides  does  not  enter  into 
the.  problem  of  this  work,  we  will  not  in  detail  analyze  this- case  and 
phenomenon  of  degeneration,  referring  to  work  [1].  Now  it  is  possible 
to  visualize  the  common  picture  of  behavior  X„,  in  the  dependence  on 
the  thickness  of  waveguide  h.  The  corresponding  trajectories  are 
schematically  represented  in  Fig.  4  and  5.  Arrows/pointers  showed  the 

direction  of  the  motion  of  root  with  increase  of  h.  Fig.  4 

'2 

corresponds  to  medium  with  em>0,  while  Fig.  5  «m<0  (cold  plasma). 

The  parts  of  trajectories  xm,  arranged/located  on  the  lower  sheet, 
are  depicted  as  dotted  line. 
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For  the  subsequent  analysis  us  it  will  be  necessary  to  estimate 
with  what  h  of  trajectory  of  zeros  the  half-lines  (1.5)  intersect.  If 
we  use  formulas  (2.6)  (2.7)  and  (2.12)  and  to  be  interested  in  upper 
half-plane  (A),  then  the  height  of  waveguide  h„,  at  which  trajectory 
m  -  zero  will  cross  half-line  (1.5),  it  will  be  located  from  the 
condition 


Po  =  P>  1, 


where  Po  —  kh0,  F(t)—  and  p0 


cos 


1  |  v  | 

7= - 7*  if  condition  p>  1  is  not 

y  p  cos  — 


satisfied. 


Thus,  in  our  approximation/approach,  which  is  determined  by 
force  (2.11)  with  p>l  by  the  inequalities 


P  >?*(?)  n^H 


-J  ,  A 
npH  m  7^=  0, 

m  =  Q  («-  <  0), 


Key:  (1).  with. 


all  a,„  simultaneously  cross  line  (1.5).  Exception/elimination  is 

'j 

only  zero  A,  when  £m>0  (see  Fig.  4).  It  means,  if  J3</30,  then  with  the 
coincidence  of  the  duct/contour  of  integral  (1.1)  with  the 
duct/contour  of  most  rapid  descent  (1.5)  in  the  case  e,»>0  do  not 
appear  deductions  at  points  x„„  of  exit  point  of  which  they  satisfy 
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condition  (2.3).  But  if  e„1<0,  then  when  /3</30  the  solution  is  present 
the  twice  infinite  sum  of  deductions.  The  waveguide,  whose 
height/altitude  satisfies  the  condition 

0<(V  (2,16) 

we  will  subsequently  call  thin. 
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rpaa.Q) 

l-  lr> 

0 

1 

30 

0.91 

45 

0,78 

60 

0,57 

60 

0.22 

90 

0 

130 

1 

135 

1.8 

150 

8.34 

170 

11,1 

180 

00 

Key:  (1).  deg. 
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In  this  case,  of  course,  it  is  assumed  that  the  conditions  for 
applicability  (2.11)  of  formula  (2.12)  are  satisfied  also.  For 
convenience  in  the  estimate  of  the  magnitude  /50  the  table  gives 
function  F(<p). 

Let  us  determine  order  of  magnitude  h0  for  the  waveguide  with  an 
infinite-conductinc  lower  wall  and  with  the  properties  of  upper  wall, 
close  to  the  properties  of  the  ionosphere  at  the  low  frequencies. 
Considering  that  the  conductivity  of  latter  3u  ~  '0'"  ■  and  f-10  kHz, 

we  will  obtain  h,»1.2  km. 

Thus,  even  in  the  case  of  SDV  the  examined  waveguide  is  in  our 
understanding  "thick". 

We  will  be  occupied  now  the  study  of  zero  equations  (2.1),  which 
satisfies  condition  (2.4).  We  will  subsequently  designate  this  zero 
through  X( 0 ) .  Let 


C-M7) 
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Taking  into  account  (2.4)  equation 
can  be  recorded  in  the  form 


(2.1)  with  sufi 


or 


(2.18) 


4  —  e,«  -h  ( e«.)Jp2  (  —  1  )2 


(2,1'J) 


For  solution  (2.19)  let  us  use  the  method  of 
approximations,  after  assuming 

(yu)o  ~ 


then  it  will  seem  that 


(*o)l  =  yr~^n  T  4r^2(  e«i)  +  0($~)  j 


Is  not  difficult  also  to  find  (K)i 


(Xt0))r-A>i  +  3Wr(eJ,  ”  (2.22) 

Fi  (em)  =S  em(em  —  l).  (2.23) 


Thus,  in  contrast  to  exit  points  Xa  and 


iciently  small  h 


successive 


.(2.20) 

(2.21) 


exit  point  X ( 0 )  is 
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located  at  finite  distance  from  the  origin  of  coordinates  and 
coincides  with  branch  point.  Knowing  the  argument  of  function 
it  is  possible  to  determine,  on  what  sheet  is  arranged/located  the 
root  X(°)  equation  (2.1).  Under  the  condition 


4-  ^  arg  (£m) 


(2,24) 


the  root  X(“)  is  arranged/located  on  the  upper  sheet.  In  turn  arg F{<=.'m) 
characterizes  position  X{°)  relative  to  half-line  (1.5). 

Page  41. 

The  value  of  arguments  jP  (=,'„)  and  F(sn)  they  make  it  possible  to 
judge  Fig.  6  and  7.  In  Fig.  6  they  are  depicted  on  plane  (em)  of 
line,  along  which  arg («<«)  =  const.  Numbers  by  the  curves  correspond  to 
the  numerical  value  of  argument.  Fig.  7  makes  the  same  sense,  are 
only  here  depicted  lines  arg  F(sm)  =  const.  Ranges  of  values  of  s,„,  with 
which  (2.24)  it  is  not  made,  are  designated  in  Fig.  6  by  shading. 
With  such  values  4,  the  root  X(°)  transcendental  equation  (2.1)  is 
arranged/located  on  the  lower  sheet  of  plane  (X).  With  the  combining 
of  the  circuit  of  integration  in  (1.1)  with  the  duct/contour  of  the 
most  rapid  descent,  which  proceeds  from  the  point  X=k2  (half-line 
(1.5)),  we  had  to  take  into  account  deduction  in  the  pole  X(°),  if 
the  latter  is  arranged/located  somewhere  on  the  upper  sheet,  with 
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exception  of  region  (a)  (Fig.  2),  or  on  the  lower  -  in  region  (a). 
Being  turned  to  Fig.  6  and  7  not  difficult  to  establish  that  if 


arg/X)  >  -J-- 

arg  Fi  (e,„)  <  -j- , 


(2.25) 


then  A(  °)  is  located  on  the  lower  sheet  in  region  (a).  But  if 

arg  F\  (e'„)  >  4-  , 

P(  M  ^  5“  (2.26) 
arg^le^X-y-w, 

then  at  X ( 0 )  is  arranged/located  on  the  upper  sheet  out  of  region 
(a).  Of  the  aforesaid  it  above  follows  that  with  fulfillment  (2.25) 
or  (2.26)  the  strain  of  duct/contour  (1.1)  must  be  accompanied  by  the 
addition  of  deduction  in  the  pole  X(®).  Uniting  (2.25)  and  (2.26),  we 
will  obtain  the  condition 

~  (2.27) 


with  fulfillment  of  which  deduction  enters  into  field  expression.  In 

Fig.  7  shading  isolated  ranges  of  values  in  which  the  deduction 

'2 

in  the  pole  X(°)  should  not  be  considered. 


On  this  we  finish  the  research  of  the  dynamics  of  the  poles  of 
transcendental  equation  (2.1).  The  obtained  results  make  it  possible 
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to  pass  to  the  conclusion/output  of  analytical  field  expressions  in 
the  thin  waveguide. 

§  3.  Field  in  the  thin  waveguide. 

We  will  consider  condition  <1.6)  as  that  made  and  let  us 
represent  field  in  the  form  of  the  following  sum: 
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Fig.  6. 
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are  encompassed  the  deductions  in  poles  As  it  was  shown  in  the 

preceding/previous  paragraph,  in  the  case  s,n>0  and  thin  waveguide 
such  deductions  must  not  be  considered.  Therefore,  if  «m>0,  then 

(3,2) 


But  if  em<0,  then  with  the  coincidence  of  initial  duct/contour  with 
the  duct/contour  of  descent  will  be  crossed  all  of  pole  >>.„•  In  this 
case  iff1  will  be  represented  in  the  form  of  the  twice  infinite  Siam 
of  normal  waves.  Let  us  consider  first  component /term/ addend  t'{J] . 
Taking  into  account  (1.6),  let  us  replace  the  Hankel  function  with 
its  asymptotic  representation  and  we  deform  the  duct/contour  of  the 
integration.  As  a  result  it  will  prove  to  be 


/ 


E 


(0 


=  I  +  P, 


4*'o  V  nr  J  Y  k2  -  V- 


F(\,  h)elXrd\% 


(3,3) 


where  F(\,h)  is  determined  by  formula  (1.2)  and  through  P  it  is 
designated  deduction  in  the  pole  X(°)  (if  the  latter  enters  into  the 
solution).  The  form  of  circuit  G  is  represented  in  Fig.  8. 

Let  us  pass  from  variable  X  to  by  the  variable  u  according  to 


the  formula 


k  =  -f-  ia3. 


(3,4) 
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Then  it  seems  that  the  coordinate  of  pole  X(a)  on  the  plane  (u)  will 
be  represented  in  the  form 

uo=  Y4-$F{tm)e  1  \  (3,5) 

Point  u=0  is  saddle  point  of  integrand  I.  In  the  vicinity  of  this 
point 

K 

j/" k\  —  X2  za  Y -A,  ue  4, 

0<argl^<~,  (3,6) 

u<0  corresponding  to  upper,  and  u>0  -  to  lower  sheet  of  plane  (X). 


Page  44. 

Let  us  assume,  furthermore, 


7kY%m—  1,  (3,7) 

where  th_  fixation  of  branch  ]/£,«  —  1  is  arbitrary.  Considering  as 
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escaping/ensuing  from  assumption  about  the  fact  that  waveguide  thin, 
let  us  represent  I  in  the  form 

/  =  A  +  /S.  (3,9) 


where 


/  J^Ldv. 

1  'Jittg  r  J  v  -f  ^  sr 

—  CO 

.3  + 

o  ^  lh+r—i —  «  /»  „,2 

^  ]^=J  77^*. 


/,= 


2iu0 


YTr 


4* 


(3.10) 

(3.11) 

(3.12) 


Integrals  (3.10)  (3.11)  are  simply  connected  with  the  expression 
for  the  function  of  weakening  Sommerf ield's  problem.  They  are 
calculated  by  standard  methods  [53  and  ere  reduced  to  the  probability  • 
integrals  of  the  complex  argument.  As  a  result  of  simple  computations 
we  come  to  the  formula 


I  in  Y  o 

wx  (sr)  =  1  -f-  'lV  sre~sr  (  tt^dv. 

'*  YTr 


(3.13) 

(3.14) 


Second  term  in  the  brackets  of  formula  (3.13)  appears  from  integral 
I j  (3.11).  By  virtue  of  (3.8)  this  value  of  the  order  of  an  error  in 
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the  computations  and  we  it  subsequently  will  not  consider. 

If  we  turn  now  to  inequalities  (2.25)— (2.27) ,  then  it  will  seem 
that  in  the  expression  for  should  be  considered  the  deduction 

under  the  condition 


Ira  ^  0, 


(3,15) 


P  taking  the  following  form: 


P  =  -  [l  —  f  (4  —  1 ) I  V  V (3. 1  ii) 

~‘kE0  ' 


Taking  into  account  (3.15)  it  is  possible  to  represent  now  sum  (3.13) 
and  (3.16)  in  the  following  final  form: 


c<i>_  k7-p* 


—  &  (iV), 


V)  (sr)  =  1+2  V  sr  e  ir  j  ev'dv. 

VTr 


(3,17) 


Page  45. 

We  agreed  to  call  value  £'!•"  the  sky  wave. 

v 

In  the  case  s„ > 0  formula  (3.17)  characterizes  complete  field  in 
the  thin  waveguide.  In  the  form  it  coincides  with  the  solution  of  the 
problem  of  Sommerfield  about  the  field  of  the  dipole  above  the 
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impedance  plane.  But  if  < 0,  then  there  should  be  taken  into  account 

still  component/term/addend  from  (3.1),  that  characterizes  the 

sum  of  deductions  in  poles  K„<  which  we  call  the  sum  of  normal  and 
quasi-normal  waves, 


(3,19) 


moreover  it  is  assumed  that  |kr|>l.  Let  us  note,  besides  the  fact 
that  the  waveguide  is  thin  and  that  with  h-»0  K« I  -  Then,  if  they 
take  place  of  inequality  (2.11),  expression  for  Um  is  simplified  and 
it  will  seem  that 

l3'20) 

This  formula  we  will  use  subsequently.  However,  value  a,„  is 
determined  by  expression  (2.12),  in  which  with  m>0  should  be  taken 
the  upper,  and  with  m=0  and  m<0  lower  sign. 

v 

By  such  shape  of  formula  (3.17)  (3.18)  and  (3.20)  they  solve  the 
problem  about  the  determination  of  field  in  the  flat/plane  thin 
waveguide.  The  following  paragraph  will  be  dedicated  to  their 
analysis. 
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S  4.  Fundamental  laws  governing  the  field  in  the  thin  waveguide. 

The  investigation  of  expression  for  E:  let  us  begin  from  case 
e»<> 0.  Then  it  seems  that 

£,  =  £i"  (4,1) 

it  suffices  to  consider  expression  (3.17).  This  expression  takes  the 
form,  characteristic  for  the  solution  of  the  problem  of  Sommerfield 
about  the  field  of  the  dipole,  arranged/located  in  the  uniform 
half-space,  limited  by  impedance  plane. 

Page  46. 

For  coinciding  the  field  in  this  problem  with  the  field,  described  by 
formula  (3.17$,  one  should  consider  that  the  properties  of 
half-space,  in  which  is  arranged/located  the  dipole,  coincide  to  the 
properties  of  the  medium,  which  limits  waveguide  on  top,  and  the 
given  surface  impedance  5  of  plane  in  Sommerf ield 's  problem  is  equal 


to 
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Only  formal  difference  from  Sommerf ield's  problem  lies  in  the  fact 
that  in  our  case  sign  Re5  is  arbitrary. 

Let  us' recall  the  fundamental  properties  of  function  w(sr), 
which,  following  the  terminology  of  Sommerf ield's  problem,  we  will 
call  the  function  of  weakening,  and  its  argument  sr  -  numerical 
distance.  We  will  first  consider  that 

|  sr  i  <  1 ,  (4,2) 

i.e.  let  us  consider  the  case  of  small  numerical  distances.  Then 

w  (sr)  =z\  -\-iV  K$r  —  2  sr.  (4,3) 

It  means,  if  the  property  of  waveguide  it  satisfies  the  condition 

•'  ImVrs>0,  (4,4) 

then  |w(sr)|  decreases  with  an  increase  in  the  distance.  However, 
with  the  fulfillment  of  reverse/inverse  inequality 

Im  V  s  <C  0  (4»5) 

'2 

|w(sr)|  increases  with  increase  of  r.  However,  concerning  argument 
w(sr),  which  we  will  call  supplementary  phase  [argw (sr)=  fJon(5r)],  the 
latter  it  increases  with  increase  of  r  under  the  condition 


Re  V  s  >  0 


(4,6) 
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and  it  decreases  with  the  reverse/inverse  inequality.  Let  us  note 
that  by  virtue  of  (3.12) 

arg  V  s  it  +  arg  (4,7) 

However,  on  the  value  arg/^iO  with  different  ones  it  is  possible 
to  judge  with  the  help  of  Fig.  7.  Being  turned  to  this  figure,  is  not 
difficult  to  conclude  that  in  the  first  quadrant  (appropriate 
•;>o)  inequality  (4.4)  has  the  place  in  region  (1)  and  in  the  part 
of  region  (6),  and  inequality  (4.6)  in  region  (4)  (5)  and  also  in 
part  (6). 


If  one  considers  that  the  phase  rate  in  waveguide  Cty  10/IM 
connected  with  phase  rate  r*  of  plane  wave  in  medium  z>h  with  the 
relationship/rat io 


(4.8) 


then  it  is  possible  to  conclude  that  in  regions  (1)  (2)  (3)  with 

\* 

small  |sr|  c$  so™  >  decreases  with  the  increase/growth,  r,  and  in 
regions  (4)  (5)  and  (6)  (partially,  since  it  is  examined  the  first 
quadrant)  c't> and  increases  with  increase  of  r. 
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From  formula  (3.17)  it  follows  also  that  if 
argVs  =  -|-*,  (4.9) 

then  <f>Ao,.  =  °  with  all,  r,  whence  follows  £«li  Honn 

Relationship/ratio  (4.9)  occurs  on  the  line,  which  divides  regions 
( 3 )  and  ( 4 ) . 

Let  us  consider  now  the  behavior  of  field  with  the  large 
numerical  distances 

|  .st  |  »  1 .  (4.10) 

In  this  case  it  is  possible  to  note  w(sr)  by  its  asymptotic 
representation 


w(sr)~ 


(2/i  -  1) !' 

2n  ( $r)n 


im  yr  s  >  0, 
'2 


21 V  Ksre~sr 


(2»  —  i)i!  . 
j  2" (sr)n  ’ 

n-~\ 


Im  V  <  0. 


(4,11) 
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Holding  only  one  member  in  expansions  (4.11),  we  will  obtain 


w  ( sr ) 


lm  Vs  >  0, 

—  b(sr);  Im  Vs  <  0, 
O  (sr)  =  1  —  M  V  re  ( sr)3/2e~ir . 


(4.12) 

(4.13) 


Most  simply  field  appears,  if  the  properties  of  upper  half-space 
are  described  by  certain  point  from  region  (1).  In  this  case  of  w(sr) 
it  behaves  just  as  in  the  case  of  the  strongly  capacitive  plane  of 
the  Earth:  |w(sr) |  monotonically  decreases  with  increase  of  r,  and 
fPAoa  approaches  the  constant,  equal  to 


Vton  (sr)  -  *  —  arg  s  =s  —  —  2  arg  F{zm).  (4,14) 

r -*■  oo 


Behavior  of  |w(sr)|,  <p«oa (sr)  and  ^<huojih  for  the  waveguide  with  the 
properties  in  question  is- represented  in  Pig.  9-11. 

Thus,  dependence  on  the  distance  of  field  £«  in  the  field 
|sr|»l  is  described  by  factor 


— r-.  (4.15) 

but  phase  rate  with  r-*°>  approaches  the  phase  rate  of  plane  wave  in 
the  unbounded  medium  with  the  properties  of  upper  half-space. 
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Let  us  now  move  on  to  the  analysis  of  field  at  the  large 
numerical  distances  in  the  case,  when  the  properties  of  upper 
half-space  are  described  by  certain  point  in  regions  2-5  (Fig.  7), 
moreover  is  examined  only  that  part  of  region  5,  in  which  em>0. 


Page  48. 

In  the  regions  Imj/S<0  indicated  and  function  w(sr)  has 
asymptotic  representation  (4.13).  We  will  be  occupied  the  study  of 
the  function  (D(ir)  from  (4.13),  which  let  us  rewrite  in  the  form 

— <  1 1 1  r  sin  *•(  .ty— < 

®(sr)  =  1  -f-p(r)  e  *,  (4,K>) 

'!>  =  arg  /"s, 
p  (r)  ==  4  v  TC I  sr|  e 

Being  turned  to  Fig.  7,  let  us  note  that  the  angle  2$  is  changed 
in  the  following  limits: 


2tc<2-k4* 

o6^acTi<0(2), 

(4,17) 

j*<2'K3  * 

owiuCTt^ip), 

(4,18) 

^  2-|>  C  -J- - 

oo/iacTt@(4), 

(4,19) 

4 * ^  4* 

oo;iacTi$;5). 

(4,20) 

Key:  (1).  region. 
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It  means  that  in  regions  (2)  and  (5) 


P^)  -0,  (4,21) 


and  in  regions  (3)  and  (4) 


p(r)  -*■  co. 
r  -*co 


(4.22) 


DOC  =  84100102 


PAGE 


Fig.  11. 

Page  49. 

Let  us  consider  case  (4.21).  Function  p(r)  has  a  maximum  at  a 
distance  in  this  case 

I  s  I  ’  (4,23) 


moreover 


(4.24) 
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It  is  not  difficult  also  to  be  convinced  of  the  fact  that  the 
equation 

~('-i0}=l  (4,25) 

has  two  roots 


|s|ri»  ~  A!J_ 

1  1  '  cos  2-4/  • 


(4,26) 


Form  p(r)  is  represented  in  Fig.  12. 

However,  as  far  as  function  #(sr)  (4.16)  is  concerned,  in  case 
(4.21)  in  question  on  the  complex  plane  ($(sr))  depending  on  r  it 
will  take  the  form  of  the  spiral,  which  emerges  from  point  +1.  with 
r-»«  point  +1  is  limit  point.  Behavior  of  $  in  region  (2)  is 
represented  in  Fig.  13.  In  region  (5)  4>(sr)  is  analogous  spiral,  only 
rotation  around  point  +1  occurs  counterclockwise.  Let  us  designate 
the  number  of  revolutions  of  spiral  around  the  origin  of  coordinates 
th  ough  m.  Then  limiting  value  (Paou  is  located  from  the  formula 


<?iun  ‘ 

r  co 


arg  s  —  2m- 1 

m=l,  2.  3,  ... 


(4,27) 
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In  turn  Mach  number  coincides  with  a  number  of  points  of  intersection 
with  the  spiral  of  the  negative  part  of  the  real  axis  of  plane  ($) . 
However,  the  latter  for  region  (2)  are  determined  by  the  following 
relationships/ratios: 

(  |  $  |  r4  sin  2'Ji  —  3<|»  +  -j-  =  (2k  -f-  1  ).ji, 

A  =  0,  1,  2  ...  (m  —  1). 
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Let  us  attempt  to  find,  with  what  properties  of  upper  half-space 
in  formula  (4.27)  m=0.  For  this  let  us  note  that  the  smallest  root  of 
equation  (4.28)  corresponds  k=0  and  is  equal  to 


(4.29) 


If  it,  seems  that 

r0>rm,  (4.30) 

then  (4.28)  will  not  have  the  solutions,  i.e.,  in  formula  (4.27)  m=0. 
Considering  (4.26)  and  (4.29),  it  is  possible  to  give  (4.30)  the 
following  form: 

x  3 

T  + 

1  -r-jj-  sin*  4. 


(4.31) 
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After  using  the  method  of  the  successive  approximations  (in  the  zero 
approximation  we  assume/set  tg2i//0=7r/8 )  ,  it  is  not  difficult  to  find 
the  values  r  which  satisfy  inequality  (4.31), 

*<*<3,39  (4,32) 


or 

-^-<arg/r(enX  1,04.  (4,33) 

The  part  of  region  (2),  in  which  is  satisfied  condition  (4.33),  we 
will  call  "weak"  part.  Here  function  w(sr)  has  depending  on  distance 
comparatively  weakly  expressed  interference  character  and,  in 
particular,  |w(sr)|  it  cannot  become  zero,  but  limiting  value 
<P/wii(sr)  is  described  by  formula  (4.27)  with  m=0.  Analogously  is 
isolated  "weak"  part,  also,  in  region  (5),  where  it  is  determined  by 
the  inequality 


-|-K>arg/r(em)>3,68. 


(4,34) 


*'ig.  14  depicts  regions  (2)  and  (5)  with  the  chosen  "weak"  parts.  It 
is  now  not  difficult  to  visualize  the  dependence  of  field  in  the 
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waveguide  on  the  distance  in  the  "weak"  and  "strong"  parts  of  regions 
( 2 )  and  ( 5 ) . 

Let  us  consider  region  (2)  first.  The  corresponding  graphs  for 
|w(sr)|,  fnon (r)  and  «ujih  are  represented  in  Fig.  15-17. 

With  the  change  in  the  properties  of  upper  half-space,  which 
corresponds  to  the  motion  of  representative  point  from  the  "weak" 
part  of  region  (2)  into  the  "strong"  (|s|=const),  the  field  endures 
the  following  changes:  a)  increases  the  amplitude  of  the  oscillations 
of  field  and  their  three-dimensional/space  period  increases;  b) 
establishment  <pwu(sr)  occurs  at  large  distances  from  the  source;  c) 
also  more  slowly  is  installed  c*  ooJ1H  to  value  <\j,. 
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C*  to** 

rt 


Fig.  17.  Fig.  18. 

Fig.  17.  Key:  (1).  (Strong).  (2).  (Weak.). 

Fig.  18.  Key:  (1).  (Strong).  (2).  (Weak.). 

Page  52. 

However,  with  r-*«  field  in  region  (2)  and  (5)  behaves  just  as  in 
region  (1),  dependence  on  the  distance  as  before  takes  form  v4.1c'. 

But  if  the  properties  of  upper  half-space  are  described  by  that 
part  of  region  (5),  which  corresponds  em>0,  then,  as  has  already 
been  indicated,  the  properties  of  fields  will  be  close  to  case 

examined  above.  The  character  of  behavior  |w(sr)|  will  remain  as 

*» 

before,  however,  concerning  <Pa«.n,  that  instead  of  formula  (4.14) 
should  be  utilized  a  relationship/ratio  1 
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(sr)  -  4-  -  -  2  arg  F  (.I) 

r  oo 


(4,05) 


in  "week"  p-ft  (5)  and 


‘PjOII  ('"^) 

r 


4-  -  —  2  arg  /=■(£«,)  -f  2/n", 
c.->  m  =  1 ,  2,  3  . . . 


(4,30) 


in  the  "strong"  ^art. 

PC  tf^TE  1 ,  It  is  assumed  that  arg  w(o)=0.  ENDFOOTNOTE. 

Dependence  %wn(sr)  in  the  region  in  question  is  represented  in  Fig. 

18,  and  nature  of  behavior  e.i. ■•«,.»(/•)  19.  Thus,  the  values  of  the  phase 
rate  and  (P»°«  in  region  (5)  in  contrast  to  (2)  with  small  |sr| 
increase  with  an  increase  in  tne  distance  from  the  source. 

.Let  us  now  move  on  to  the  examination  of  field  in  the  case,  when 
the  properties  of  upper  half-space  are  described  by  certain  point 
from  region  (3)  or  (4)  Fig.  7.  The  behavior  of  field  with  small 
numerical  distances  was  already  considered  above.  We  will  therefore 
consider  that  |sr|>>l.  Let  us  note  that  in  these  regions  an<^ 

at  the  large  numerical  distances  from  the  source  in  expression  (4.13) 
should  be  considered  only  second  component/term/addend.  Then  it  seems 
that  (3.17)  it  will  be  possible  to  reduce  to  the  simpler  form 
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£. 


k  iP o  _  _  r - 


'lV  ~s  — 


ik,r  —  if  :  <  — 


V  r 


(4,37) 


If  we  compare  the  dependence  of  field  on  the  distance,  described  by 
expression  (4.37),  with  the  previously  cases  (formula  (4.15) 
examined),  then  it  is  not  difficult  to  conclude  that  in  regions  (3) 
and  (4)  occur  following  properties  £;. 
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Fig.  19.  Key:  (1).  (Weak.).  (2).  {Strong). 

Page  53. 

i 

a)  divergence  is  described  by  factor  ~7f  >  which  is 
characteristic  for  the  waveguide  problems.  However,  in  the  regions 
(1)  (2)  and  (5)  the  corresponding  factor  took  the  form  -L  .  This 
dependence  takes  usually  place  in  the  problem  from  the  field  of  point 
source,  located  near  interface  in  the  two-layered  (but  not 
three-layered)  medium. 

b)  exponential  factor  takes  now  the  form 

r**'-4'.  (4,38) 

This  means  that  the  limiting  value  of  phase  rates  now  no  longer 
coincides  with  the  phase  rate  of  plane  wave  in  upper  half-space  c\t, 
(which  is  also  characteristic  for  Sommerf ield’s  problem),  and  it  is 
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determined  by  the  expression 

U‘|.  »o."iJii|i<4  —  _  |  i  i  ->iM  ■ 

1  “  Re  ft. 

It  means,  in  region  (3)  and  in  region  (4)  c . . The 

behavior  of  the  phase  rates  in  the  dependence  on  the  distance  is 
represented  in  Fig.  20  {region  (3))  and  21  (region  (4)).  Furthermore, 
if  in  all  previously  cases  in  question  the  exponential  character  of 
the  decrease  of  field  with  the  distance  coincided  with  the  decrease 
of  the  plane  wave,  which  is  propagated  in  the  medium  with  the 
properties  of  upper  half-space,  then  in  this  situation  the 
corresponding  factor  takes  the  form 


g- |lm*,-|»|  |  r 


(4,38a) 


This  means  that  the  field  in  the  waveguide  decreases,  generally 
speaking,  it  is  slower  than  field  in  the  upper  half-space.  This 
difference  most  strongly  is  expressed  on  the  line,  which  divides 
regions  (3)  and  (4).  Thus,  it  proves  to  be  that  in  the  case  of  media 
with  properties  (3)  or  (4)  is  observed  more  clearly  the  expressed 

'X 

channeling  of  field.  This  is  connected  with  the  less  effective 
penetration  of  the  latter  into  the  upper  half-space. 


Let  us  now  move  on  to  the  examination  of  the  case,  \vhen  upper 
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half-space  represents  conducting  medium  with  the  negative  dielectric 
constant.  As  has  already  been  indicated,  cold  plasma  possesses  such 
properties. 
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Fig.  20.  Fig.  21. 


Page  54. 


For  describing  the  field  in  the  waveguides  with  similar  properties  it 
is  necessary  in  expression  E,  (3.1)  to  consider  also  component  £t21  of 
(3.18)  —  (3.20) .  If  we  use  formula  (2.12),  then  will  be  represented 
in  the  form 


>c  f  |  v  |  —  r 

Cz  ~  2«o  r 


e  R  W,  (r). 


(4,39) 


TVt(r) 


—  |/" _L_  |[  |  v  j  -r  t«|J  *  +  [i  (mi t  —  u)  -f  |  v !  j  12  > 

*). ' 


-  (m*-2u)  — I 


(4,40) 


where  u  and  v  are  assigned  by  formula  (2.7),  moreover  by  virtue  of 

'i 

(2.8)  u^7.72.  As  it  follows  from  formulas  (4.39),  all  waves  entering 
expression  E?\  move  with  identical  phase  rate  equal  to 
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where  c  -  speed  of  light  in  the  vacuum,  The  effective  depth  of 
penetration  of  such  waves  into  waveguide  d?  can  be  found  from  the 
formula 


(4,42) 


If  leml~c,(°).  then,  as  is  evident,  for  example,  from  Fig.  3  u~ir/2 
and  As  we  will  see  subsequently,  in  this  case 

component/ term/ addend  £l*’  very  little  in  comparison  with  £[l).  But  if 
|41»1,  then 


it  proves  to  be  that 


arctg-^-L 


0 

e0 

|2  '  > 


V : 


i  >  Kn 


(4,-13) 


i.e.  the  depth  of  penetration  can  considerably  exceed  the 
thickness  of  waveguide  h. 
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Let  us  now  move  on  to  the  computation  of  component/term/addend 
£”>  from  (3.1).  Let  us  preliminarily  note  that  value  \p  from  (4.16)  in 
the  appropriate  regions  of  the  second  quadrant  of  plane  is 

determined  by  inequalities1 

— <cf’o6JiacTb  (5), 

0  ^  •}<  ^  o^acTb  (p 

4  ,f)  (4,44) 

—  <  +  <  —  oo^acTb  (7). 


Key:  (1).  region. 

FOOTNOTE  1 .  Here  from  both  parts,  of  the  iue».  quality  unessential  for  us 
value  2jt  is  subtracted.  ENDFOOTNOTE. 

Page  .5.5. 

It  means,  for  determination  E[l)  in  the  part  of  region  (5),  which 
lies  in  the  second  quadrant,  should  be  used  formula  (4.37),  and  in 
regions  (6)  and  (7)  -  by  formula  (3.17).  Let  us  note,  besides  the 
fact  that  under  the  condition.? 
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makes  sense  to  examine  only  regions  (6)  and  (7).  However,  in  the 
latter  dependence  Et  on  the  distance  with  the  the  large  |sr|  is 
determined  by  factor  (4.15) 

elk,r 


Consequently,  for  the  part  of  the  field,  described  by 

component /term/ addend  £*",  effective  depth  of  penetration  is  equal  to 


4l> - L__. 

k  I  in  V i'm 


■  (4,46) 


and  in  the  cac-  cf  small  ones  in  question  !  £m !  considerably  exceeds 
</<*>..  Thus,  case  (4.45)  is  reduced  to  study  £*-l)  in  regions  (6)  and 
( 7 ) *  Taking  into  account  (4.44),  it  is  not  difficult  to  conclude  that 
here  lm/s>0,  and  this  means  that  the  field  behaves  qualitatively  just 
as  in  region  (1).  Difference  consists  only  of  the  fact  that  here  in 
view  of  Re/s>0  c<\ c  from  the  side  of  smaller  (but  not  high) 

r  —  ~ 

values.  The  oscillations  of  phase  rate  also  are  not  observed. 


Let  us  now  move  on  to  the  examination  of  the  more  interesting 


case 
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K!»>. 

Hosm<0, 


when1  can  prove  to  be  that 


or 


45)  -4“ 

d?  »  4“. 


(4,47) 


(4.48) 

(4.49) 


FOOTNOTE  1 .  In  the  "asymptotic"  thin  waveguides  h-»0,  o  and 
cannot  be  thus.  ENDFOOTNOTE. 

Page  56. 

If  it  seems  that  inequality  (4.49)  occurs,  then  at  the  distances 
from  source  r>dlcl)E.^£ ‘2|.  However,  value  E?]  is  defined  as  the  sum 
of  waves  with  the  cylindrical  divergence,  i.e.,  by  "waveguide"  type 
waves.  Thus,  in  contrast  to  previously  case  (4.45)  examined  here  is 
observed  the  more  clearly  expressed  channeling  of  energy.  This 
channeling  was  observed  when  em>0  only  in  regions  (3)  or  (4).  In 
the  case  (4.46)  components/terms/addends  w”  and  E?'  are  equally 
essential.  Since  phase  wave  velocity,  which  correspond  to  these 
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components/terns/addends ,  are  different,  complete  field  will  carry 
complicated  interference  character.  Especially  complicated  form  will 
have  a  field  in  region  (5),  where  as  a  result  of  the  presence  of 
deduction  P  (3.16)  interaction  of  three  waves  with  different  phase 
rates  will  be  observed.  However,  if  together  with  (4.48)  is  satisfied 
condition  then  sufficiently  far  from  the  source 

~.(2) 

interference  will  not  be  observed,  since  component/term/ addend 
has,  as  has  already  been  indicated,  cylindrical  divergence,  but 
divergence  E(zl)  is  determined  by  factor  l/r2.  Thus,  in  this  case  at 
large  distances  from  the  source  the  "waveguide"  character  of 
propagation  as  before  occurs. 
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SUPPLEMENT  AND  REVISION  TO  THE  ARTICLE  "On  THE  TROPOSPHERIC 
REFRACTION  OF  RADIO  WAVES"  [1]. 

G.  I.  Makarov,  N.  P.  Tikhomirov. 

In  work  [1]  the  problem  of  diffraction  radiowave  propagation 
around  the  Earth  upon  consideration  of  the  atmosphere  heterogeneous 
in  the  height/altitude  was  examined  and  the  approximate  analytical 
expressions  for  the  zero  eigenvalue  t0  were  obtained. 
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In  particular,  with  the  sufficiently  small  surface  impedance  of 
formula  fort,  they  take  the  following  form  1 : 


with  the  linear  profile/airfoil  of  the  atmospheric  heterogeneity 


'/> 


‘o  -*0 


JJi _ 3_  .. 

—  *l» 


(1) 


with  the  exponential  profile/airfoil  at  the  high  l) 

frequencies 


_i_  r.u  \:i _ 

t  __  I  I.  ,  3  3.5  \  oj  10 

*0— A)-  —  +  ^  +  - T, - |is,  (2) 


t"  1  2/;;  ‘  i  + «  3// 


.  with  the  exponential  profile/airfoil  at  the  low  (*i£l) 
frequencies 


2/“ 


2  +  (A3 


1  + 


(3) 


FOOTNOTE  1 .  Here  we  do  not  dwell  on  the  designations,  they  completely 

'j 

coincide  to  designations  in  [1].  ENDFOOTNOTE. 


First  term  in  the  right  side  in  (1)  corresponds  to  the  concept 
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of  an  equivalent  radius,  second  term  appears  due  to  the  modified 
impedance,  and  the  third  -  due  to  Schwarz's  derivative  in  the 
equation  for  the  eigenfunctions  of  problem.  First  two  terms  in  (2) 
make  the  same,  sense,  that  also  in  formula  (1),  third  composed  here 
appears  as  a  result  of  a  difference  in  the  exponential 
profile/airfoil  from  the  linear.  In  formula  (3)  first  term 
corresponds  to  homogeneous  atmosphere,  the  second  is  caused  by  the 
effect  of  atmospheric  heterogeneity  on  surface  impedance,  and  the 
third  considers  tropospheric  refraction.  For  the  evaluation  of  the 
role  of  each  component/term/ addend  in  (1)  (2)  in  [1]  phase  rate  and 
attenuation  of  zero  normal  wave  were  calculated;  however,  during  the 
computations  was  allowed  the  error,  as  a  result  of  which  the  numbers 
given  in  Table  1.2  in  [1]  do  not  correspond  to  formulas  (1)  (2), 
i.e.,  to  formulas  (49)  (55)  on  numbering  [1]. 
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Table  1. 


/.  K.’l 

0) 

I, 

<<n  _ 

a  1  6 

II 

a 

6 

o 

0,020 

0.930540  j  0,080275 

0.930299 

3. 431.  10  - 

(  >38-  nr  - 

•>.0 

o,i)()(>:> 

0.005744  '  0,90573 1 

0  005731 

7.304- 10  J 

7.381  -10'  1 

lO.o 

0.0040 

0.097314  i  0.007300 

0  007300 

4,(184-  10"* 

4(155- i<r:i 

50,0 

0,0014 

0,000080  '  0,000070 
'  1 

0.009079 

1,59:5-  10 

1.504  10"' 

I  ’  i  I 


Note:  a)  without  taking  into  account  second  and  third  terms  in  (1), 
b)  taking  into  account  of  second  term,  c)  taking  into  account  second 
and  third  terms.  The  constant  c  during  the  computations  was 
assumed/set  equally  3*10*  to  m/s. 

Key :  ( 1 ) .  kHz . 

Page  203. 

Table  .1,  2  with  the  corrected  values  of  rate  and  attenuation  of  zero 
normal  wave  when  &„■  =  0.  are  given  below.  Furthermore,  for  the 
evaluation  of  the  role  of  third  component/term/ addend  in  (3)  here  is 
given  Table  3,  in  which  for  the  series/row  of  frequencies  are  given 
phase  rate  and  attenuation  of  zero  mode.  Table  1  corresponds  to 
linear  profile/airfoil,  moreover  in  contrast  to  [1],  into  it  are 
introduced  also  the  values  of  the  phase  rate  upon  consideration  of 
Schwarz's  derivative. 
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From 'Table  1  it  is  evident  that  with  the  linear  profile/airfoil, 
which  corresponds  to  the  lower  layers  of  the  earth's  atmosphere,  at 
frequencies  of  higher  than  5  kHz  the  effect  of  the  modified  impedance 
on  the  phase  rate  affects  only  in  the  fifth  significant  place,  and 
attenuation  -  in  the  third.  The  effect  of  Schwarz's  derivative  on 
these  frequencies  is  generally  is  negligibly  small. 

A  difference  in  the  exponential  profile/airfoil  from  the  linear 
at  frequencies  of  higher  than  10  MHz,  as  can  be  seen  from  Table  2, 
virtually  does  not  affect  the  phase  rate  of  zero  normal  wave,  but  in 
the  attenuation  it  is  developed  only  in  the  third  significant  place. 
In  Table  2  case  b,  which  corresponds  to  the  account  of  second  term  in 
(2)  is  not  reflected;  its  contribution  to  the  eigenvalue  at  the 
frequencies  in  question  can  be  disregarded/neglected.  The  role  of 
tropospheric  refraction  at  the  low  (m£1)  frequencies  is  visible  of 
given  below  Table  3. 

The  effect  of  the  modified  apparent  resistance  on  the  phase  rate 
and  the  attenuation  of  zero  mode  at  frequencies  is  above  10  kHz,  as 
are  shown  calculations,  the  second  order  of  smallness  in  comparison 
with  the  effect  of  tropospheric  refraction;  therefore  in  table  3  case 
b,  which  corresponds  to  the  account  of  second  term  in  (3)  is  not 


reflected. 
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Table  2. 


/.  Mtn 

_ mJ 

'S 

11)1 

\k 

.0 

»*  K 

a 

a 

U 

! 

10,0 

0.003 

o  999973 

0.9‘)997:i 

4.664-10—'’ 

4,704- 10— 5 

20 1) 

0,039 

0,999983 

0,999980 

2.938- 10  J 

2,954-  10~5 

50.0 

0.0?  1 

0  999991 

0,999991 

1,595- 10-5 

1,500- 10~5 

Note: 


a) 


for  linear 


prof ile/airfoil. 


c)  for  the  exponential. 


Key :  ( 1 ) .  MHz . 


Table  3. 


/.  /..!( 

\* 

J 

II 

“  i 

10,0- 

0,991)671; 

0.990040 

5,770- ID — ‘ 

5,718-10-' 

30,0 

2,71 

0,998.199 

0  998310 

2,777  - 10 — J 

2,083  10"' 

50,0 

1.9! 

0,998801 

0  998727 

1,975- 10— ‘ 

1,804- 10—' 

100,0 

1,21 

0,999282 

0,999053 

1,244-10-' 

1,188-10-' 

Note:  a)  for  homogeneous  medium,  c)  taking  into  account  third 
component/ term/addend  in  (3). 


Key:  (1).  kHz. 
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In  conclusion  the  authors  express  appreciation  to  L.  N.  Lutchenko, 
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that  indicated  the  nonconformity  of  those  given  in  [li  Table  1,  2  to 
Formulas  (49)  (55). 
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